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Design of a Surface Deformation Measuring Instrument for 
the Surface Tension Driven Convection Experiment (STDCE-2) 


1.0 Introduction 

This Final Technical Report covers the work accomplished (under NAG3-1300) from 1 
October 1991 to 1 October 1993. The grant is a direct result of Dr. H. Philip Stahl's 
(of Rose-Hulman Institute of Technology) participation in the NASA/ASEE Summer 
Faculty Fellowship Program at NASA Lewis Research Center sponsored by Case 
Western Reserve University and the Ohio Aerospace Institute. 

The Surface Tension Driven Convection Experiment (STDCE) is a fundamental fluid 
physics experiment designed to provide quantitative data on the thermocapillary flow 
of fluid under the influence of an increased localized surface temperature. STDCE 
flew on the Space Shuttle Columbia in the First United States Microgravity Laboratory 
(USML-1) in June 1992. The second flight of this experiment (STDCE-2) is scheduled 
for 1995. The specific science objectives of STDCE-2 are to determine the extent and 
nature of thermocpillary flows, the effect of heating mode and level, the effect of the 
liquid free-surface shape, and the onset conditions for and nature of oscillatory flows. 
In order to satisfy one of these objectives, an instrument for measuring the shape of 
an air/oil free-surface must be developed. 


2.0 Objectives 

The purpose of the first year of this project was to conduct a feasibility study to 
develop one or more non-contact (i.e. optical) techniques for measuring the surface 
deformation of an air-oil free-surface under the influence of a localized thermal load 
(i.e. C02 laser beam or electric heater) which could be integrated into the anticipated 
reflight of the Surface Tension Driven Convection Experiment (STDCE-2). While 
initially these measurements were to be performed in a laboratory environment, the 
instrument had to be capable of operating in a Space Transport System (STS) 
environment for the United States Microgravity Laboratory flights USML-2. 

The primary purpose of the second year of this project was to design, build, calibrate, 
and fully characterize the instrument(s) prototyped during the first year. A secondary 
purpose was to actively participate in the design of an engineering model and flight 
system for STDCE-2/USML-2. 

Additionally, both years had the goal of developing algorithms and software to 
automatically analyze the output data provided by the surface deformation measuring 
instrument. 



3.0 Accomplishments 


A prototype instrument was developed, designed, built, calibrated, and fully 
characterized for measuring the surface deformation of an air-oil free-surface under 
the influence of a C02 laser beam. The prototype instrument utilizes a modified two- 
channel Ronchi test technique where a collimated beam of light is reflected off of the 
oil free-surface and passed through a sinusoidal amplitude grating. The transmitted 
irradiance pattern provides quantifyable data about the surface's slope as a function of 
position. These results are described in great detail in the following publications: 

Meyers, William S., Design and Calibration of a Two Channel Ronchi System to 

Contour a Dynamic Fluid Surface . Master's Thesis, Rose-Hulman Institute of 

Technology, 1993. 

Meyers, Wiliam S., and H. Philip Stahl, "Sensitivity of Two-Channel Ronchi 

Test to Grating Misalignment," SPIE Vol. 1994 . 1993. 

Meyers, Wiliam S., and H. Philip Stahl, "Contouring of a Free Oil Surface," 

SPIE Vol. 1755 . 1992. 

Meyers, Wiliam S., Brent C. Bergner, Ronald D. White, David J. Huff, and H. 

Philip Stahl, "Contouring of a Free Oil Surface," SPIE Vol. 1779 . 1992. 


Algorithms to analyze the output slope data provided by the surface deformation 
measuring instrument were developed and reported in: 

Fischer, David J., John T. O'Bryan, Robert Lopes, and H. Philip Stahl, "A 
Vector Formulation for Interferogram Surface Fitting," Applied Optics . Vol. 32, 

No. 25, pp.4738, 1 September 1993 

Fischer, David J., and H. Philip Stahl, "Mechanism for Surface Fitting of 
Interferometric Slope Data," SPIE Vol. 2003 . 1993. 

Fischer, David J., and H. Philip Stahl, "A Vector Formulation for Ronchi Shear 
Surface Fitting," SPIE Vol. 1755 . 1992. 

Fischer, David J., John T. O'Bryan, and H. Philip Stahl, "A Vector Formulation 
for Interferogram Surface Fitting," SPIE Vol. 1779 . 1992. 

The developed algorithms were implemented in a software program called THIN which 
has been delivered to NASA Lewis Research Center (LeRC). This program is being 
extended by NASA LeRC to automatically analyze Ronchi patterns using image 
processing techniques. 



A modular flight system was designed for the anticipated reflight of the Surface 
Tension Driven Convection Experiment (STDCE-2) on the United States Microgravity 
Laboratory (USML-2). It was designed to satisfy the following science requirements: 

• observe the entire surface for three different sized chambers 

(12 mm, 20 mm, and 30 mm). 

• measure slopes as small as 5 pm/mm. 

• measure slopes as great as 30 pm/mm without vignetting. 

• spatially sample the surface at least 10 times per diameter. 

• reconstruct the surface shape along a given diameter. 

• and, indicate when the oil is filled to its flat surface position. 

The design of this instrument is published in: 

Stahl, H. Philip, "Conceptual Design of a Surface Deformation Measuring 
Instrument for the Surface Tension Driven Convection Experiment (STDCE-2)," 
SPIE Vol. 2005 . 1993. 


4.0 Student Involvement 


Aside from the technical objectives, a major goal of this project was to provide a 
number of students with their first exposure to a ’real-world' research project. This 
was accomplished and the Principle Investigator is exceedingly proud of each 
student's achievements. 

William S. Meyers developed and calibrated the prototype instrument for his Master's 
of Science in Applied Optics with assistance from Brent Bergner, Leif Sorensen, Ron 
White, and Dave Huff. Bill is now pursuing a Ph.D. in Bio-Physics at the University of 
Indiana/Purdue University in Indianapolis. He is performing research on optical 
tomographic imaging of the human body. 

David J. Fischer developed and encoded the Ronchi pattern slope analysis algorithms 
with assistance from John O'Bryan and Dr. Robert Lopez of the Mathematics 
Department. Dave was the recipient of a full fellowship and is pursuing a Ph D. in 
Optics at the University of Rochester 

The software package THIN was developed through the efforts of: Kent Flint, Brad 
Freriks, Kurt Louis, Tony McAllister, Steven Reid, Brad Rodeffer, John Snider, Jason 
Snyder, Kevin Stultz, and Dr. Dale Oexmann of the Computer Science Department. 

Kevin Stultz is pursuing an M.S. in Computer Science at the University of Alabama in 
Huntsville. He is performing research on image processing methods to analyze 
Ronchi patterns. His research is a direct extension of the work begun under this 
grant. 

The participation of each student is summarized in the attached table. 

Finally, two students participated in the NASA/ASEE Summer Internship Program at 
NASA LeRC sponsored by OAI and CWRU - Bill Meyers (Summer 92 and 93) and 
Kevin Stultz (Summer 93). And, Dave Fischer participated in the NASA/ASEE 
Summer Internship Program at NASA Langley in the Summer of 1992. 






5.0 Master's Thesis 


Meyers, William S., Design and Calibration of a Two Channel Ronchi 
System to Contour a Dynamic Fluid Surface . Master's Thesis, Rose- 
Hulman Institute of Technology, 1993. 


6.0 Publications 

Fischer, David J., John T. O'Bryan, Robert Lopes, and H. Philip Stahl, "A 
Vector Formulation for Interferogram Surface Fitting," Applied Optics . 
Vol.32, No. 25, pp.4738, 1 September 1993 

Fischer, David J., and H. Philip Stahl, "Mechanism for Surface Fitting of 
Interferometric Slope Data," SPIE Vol. 2003 . 1993. 

Fischer, David J., and H. Philip Stahl, "A Vector Formulation for Ronchi 
Shear Surface Fitting," SPIE Vol. 1755 . 1992. 

Fischer, David J., John T. O'Bryan, and H. Philip Stahl, "A Vector 
Formulation for Interferogram Surface Fitting," SPIE Vol. 1779 . 1992. 

Meyers, Wiliam S., and H. Philip Stahl, "Sensitivity of Two-Channel 
Ronchi Test to Grating Misalignment," SPIE Vol. 1994 . 1993. 

Meyers, Wiliam S., and H. Philip Stahl, "Contouring of a Free Oil 
Surface," SPIE Vol. 1755 . 1992. 

Meyers, Wiliam S., Brent C. Bergner, Ronald D. White, David J. Huff, 
and H. Philip Stahl, "Contouring of a Free Oil Surface," SPIE Vol. 1779 . 
1992. 

Stahl, H. Philip, "Conceptual Design of a Surface Deformation Measuring 
Instrument for the Surface Tension Driven Convection Experiment 
(STDCE-2)," SPIE Vol. 2005 . 1993. 


7.0 Presentations 

Stahl, H. Philip, "Measuring Deformations in Space," University of 
Arizona Optical Sciences Center Colloquium, Tucson, AX 1993. 
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ABSTRACT 
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tracer particles (which have been uniformly mixed into the fluid) as they move through a sheet or beam of light. 
But, while STDCE had no means for measuring surface deformation, STDCE-2 will (Figure 2). 

A deformation measuring capability has been added to STDCE-2 because current research indicates that both 
steady-state and oscillatory flow states produce deformations of the air/oil free surface. Additionally it is 
hypothesized that there is a coupling between oscillatory flow type and deformation shape, velocity flow, and 
surface temperature field. Depending upon the experimental conditions, steady state flow produces two 
char acteristic steady state surface deformations: a ’Pimple’ (bump), and a Dimple (ho e). sci atory ow 
occurs with either pimples or dimples and produces both standing wave and rotary surface deformations. 

3.0 EVALUATION OF ALTERNATIVES 

The purpose of this project is to design, build, calibrate, and fully characterize ooeor more non-contact (i.e. 
optical) prototype instruments for measuring the surface deformation of an air-oil free-surface under the 
influence of a localized thermal load (i.e. C02 laser beam or electric heater) which can b^tegrat^ mto the 
anticipated re-flight of the Surface Tension Driven Convection Experiment (STDCE-2). While initially these 
measurements will be performed in a laboratory environment, the instrument must be capable of opera ting ; in a 
Space Transport System (STS) environment for the United States Microgravity Laboratory flights USML-2. 

Because the surface is a fluid, it cannot be contoured mechanically. It must be contoured optically. There are 
two wavs to measure a two-dimensional reflective surface: interferometncally or geometrically. For most 
optical applications, interferometers are generally preferred to geometrical tests. They are more accurate and 
aeositive^And, because they measure a surface’s shape while most geometrical tests measure a surface s slope, 
they are easier to interpret. Unfortunately, because the anticipated surface deformations are on the order of 1 
pm, conventional visible interferometry is not practical. The fringes density will be too great for Memento 
resolve. One solution to this problem is to use an infrared interferometer, but such a technique might interf 
with the I'y-giiral heating process. Therefore, a geometrical technique is required. 

Because geometrical techniques measure the surface’s slope instead of its height, they have specific advantages 
and disadvantages. They are insensitive to tilting of the free-surface cause by G-Jitter. But they require two 
orthogonal data sets in order to reconstiuct its shape. These two data seta can be either x-slope and y-slope or 
i-slope and 6-slope. Whether two data sets are truly required depends upon whether the surface is rotations > 
symmetric. If the surface is rotationally symmetric, then it may be possible to satisfy the mission s science re- 
quirements with just x-slope or r-slope. But, if the surface is non-rotationally symmetric, then it ,s nec^ry to 
acquire two orthogonal data sets. Unfortunately, the surface deformations observed in the aboratory are n 
rotationally symmetric. But, a unique feature of these deformations which may eventually lead to a simple 
solution is that their outer edge is at a constant height. 

Many classical geometric test techniques (Foucault knife-edge test, wire test, Hartmann screen test, Ronchi ^test, 
etc.) were considered. Of these, three which provide data over the entire aperture were selected for further 
evaluation: Ronchi grating, Hartmann test, and Projected grid. 

In the projected grid approach, a linear or crossed grid pattern is placed directly in front of an extended light 
source (illuminated diffuser screen). This pattern is then imaged onto the reflective sample surface. The 
reflected wavefront is re-imaged by another optical system onto the video camera. The projected gnd approac 
was eli minated because it lacked sensitivity, dynamic range, and flexible spatial sampling. 

In the single-channel Ronchi grating test, a collimated beam of light is reflected off of the sample surfaw 
(Figure 3) The reflected wavefront is focused through a grating to form a Ronchigram which ^imaged ontoa 
video camera. The Ronchi approach has many advantages. By selection of the grating spat.al ^W’*** 5 
variable sensitivity and dynamic range. By selection of the grating location, it has adjustable spat.al sampling. 
Unfortunately it does not provide simultaneous orthogonal slope data sets. 



The Hart mann screen test can provide both x- and y-dope information simultaneously. But, it has limited 
gp fltjal resolution. And, for this application, it is not practical because the screen would obscure the upper 
surface flow visualization function. An alternative is to use crossed gratings in a Ronchi test. However, the 
complex patterns produced are difficult to analyze. 

Another way to obtain simultaneous x- and y-slope data is with a two c h a nn el Ronchi approach (Figure 4). 
With the two-channel Ronchi approach, a collimated beam of light (laser or white) is projected onto the 
reflective sample surface. The reflected beam is collected, split into separate x- and y-slope measurement 
channels, passed through vertical and horizontal gratings, and projected onto an observation screen for viewing 
by a video camera. A two channel Ronchi instrument was demonstrated which simultaneously acquires x- and 
y-slope data for an input wavefront \ Unfortunately, because of packaging, cost, and schedule issues, a two 
channel system could not be implemented. 

Based upon managerial issues (packaging, cost, and schedule) and science issues (sensitivity, dynamic range, 
and spatial sampling) a single channel Ronchi system was selected. The managerial constraints and science 
requirements were then iterated, resulting in Table 1. 


'j *'*?*'. Table 1* . STDCE-2 Surface Deformation Measuring Instrument 

; Science Requirements 

Observe the entire surface for three different sized chambers: 12 

mm, 20 mm, and 30 mm,. 

- / Heasure.? slope as small as 5 f/m/mm. 

' Measure * slope as great as 30 jim/zm without vignetting. 

lillM" Spatially the surface at ' least 10 times per diameter . 

# Temporal ly sample the surface at least 30 times per second. 

Reconstruct the surface shape at least along a diameter. 

Indicate when the oil is filled to its flat surface position. 


4.0 THE RONCHI TEST 

The Ronchi test is performed by observing the shadow pattern produced when a focused beam of light is passed 
through a periodic structure, such as a line grating, which is located either at focus or away from focus a 
distance z (Figure 5). If the beam is ideal, the shadow pattern looks like the grating. But, if the beam is 
deformed, then the pattern is deformed. A Ronchigram is obtained when the shadow pattern is recorded at an 
image of the surface under test, such that it is superimposed (’painted’) upon the image. Information about the 
surfece is obtained by correlating shadow lines with physical locations on the surface. The distance between 
adjacent shadows indicates how much the surface slope changes between those two locations. The amount of 
slope required to go from one shadow to its neighbor is called the equivalent wavelength. Alternatively, the 
Ronchi test can be though of as a sheared interferogram produced by overlapping diffraction orders. 


The performance of the Ronchi system with an amplitude sine grating is completely specified by its equivalent 
wavelength: 

1 « — 

m 2(F/#) 

Since the F/# of the measurement beam is fixed by the optical system, everything depends upon grating spacing, 
d. This dependence is extremely important. It allows for the measurement sensitivity to be varied on orbit by 
jnamting into the beam gratings with different spatial frequencies (Figure 6). If the surface deformation is 
large, insert a coarse grating. If it is small, use a fine grating. 

Accuracy defines die uncertainty of a given measurement. For the current sampling system, the measurement 
accuracy is approximately X^/8. A measurement system can report numbers to an arbitrary precision, but they 
may not be accurate or repeatable. The minim um measurable value of a given system is defined to be twice its 
meas ur emen t accuracy. For the current sampling system the minimum measurable 6lope is approximately 
X^/4. The maximum amount of full chamber slope which the system can measure is determined by how many 
shadow lines (fringes) the video camera can clearly resolve in the Ronchigrem. For the current sampling 
system, the Ronchigram will be limited to approximately 16 lines pairs for a maximum slope of 16 X*,. 


Table 2. Performance Summary for F/2.25 System 


1 Grating 

I Frequency 
1 Pp/mm] 

Grating 
line Space 
[mm/lp] 

V 

[/im/mm] 

Accuracy 

[pm/mm] 

Minimum 

Slope 

[jtm/mm] 

Maximum 

Slope 

pnn/mm] 

I 0.5 

2.0 

444.44 

55.56 

uni 

7111 

1.0 

1.0 

222.22 

27.78 

55.56 

3556 

2.0 

0.5 

111.11 

13.89 

27.78 

1778 

2.5 

0.4 

88.89 

11.11 

22.22 

1422 

5.0 

0.2 

44.44 

5.56 

11.11 

711 

10.0 

0.1 

22.22 

2.78 

5.56 

356 


5.0 FLIGHT SYSTEM CONCEPTUAL DESIGN 

The flight free-snrface deformation measuring instrument was designed in six modules (Figure 7): light source, 
projection/imaging, polarization, pupil relay, measurement, and camera. Each module overcomes various 
engineering and/or packaging challenges while allowing the total system to meet the science requirements. Each 
ip/winlp is designed to be independently assembled and bench aligned before integration into the final system. 
And, each module is designed to prevent unnecessary exposure of the crew to light radiation. 

5.1 The Oil Surface 

The object under test is a free surface of 2 cSt silicon oil. It has a refractive index of 1.39. And, an irradiance 
reflection coefficient of 2.7%. Because of this low reflection coefficient, the source must have at least 10 mW 
of power for the video camera to observe the Ronchi pattern. Also, ghost reflections from the optical compo- 
nents can obscure the oil surface reflection. 



































Because the free-oil surface is a specular reflector, it can be treated as a plane mirror. To implement the 
Ronchi test, a collimated ill uminati on beam is required. Two illumination geometries were considered: Off- 
axis, and on-axis. The on-axis geometry was selected for several reasons. First, it gets the lens closer to the 
fluid surface, thus it can accept higher slope errors. Second, it has no off-axis projection distortion or cosine 
scaling. And third, it is easy to package. The on-axis geometry requires the beam to travel in two directions 
through the projection/imaging module. This can cause ghost reflection problems. 

5.2 Light Source Module 

The source modul e provides a collimated beam of light which uniformly illuminates the object under test. This 
is accomplished with a laser diode and a Galilean beam expanding telescope (Figure 8). A laser diode was 
•elected because it is a very-bright source capable of surviving the launch. (White-light and LED sources were 
also considered.) Coherence is not a factor in selecting a light source. A Galilean beam expander was selected 
because it is compact and does not have a spatial filter which could become misaligned during launch. The 
beam diameter is defined by the projection/imaging module’s entrance pupil to be 22 mm in diameter. The 
hum is magnified (30X) and truncated to balance uniformity and power. 

5.3 Projection/Imaging Module 

The projection/imaging module has two functions: illumination and imaging. Its illumination function 
ma gnifies or de-magnifies the light beam to illuminate the oil surface for each of the three chamber sizes (12 
mm, 20 mm, 30 mm). Its imaging function has two requirements: to form an image of the different sized 
surfaces under test at a fixed sized pupil location, and to pass a reflected wavefront with at least 30 /im/mrn of 
dope without vignetting. The illumination and imaging functions determine the magnification properties of this 
m~tni* The vi gnetting requirement determines its F/#. Packaging issues determines its clear aperture. 

lie illumination function is accomplished with two AFOCAL lens pairs (Figure 9), an outer pair and an inner 
pair. The outer pair forms a Keplerian telescope with unit magnification. The inner pair forms a Galilean 
telescope which is positioned about the outer pair's internal focus and functions as a pseudo field lens. This 
lens provides magnification/de-magnification when it is flipped. To illuminate the 20 mm chamber, remove the 
flip leas s uch that the 22 mm source beam is relayed unchanged. To illuminate the larger or smaller chambers, 
insert the flip lens either forward or backward to magnify or de-magnify the beam (Table 3). The final 
illumination hretp size depends upon the flip lens magnification factor. The flip lens is not symmetric about the 
focal point and docs not rotate about the focal point. 


Table 3 Image Magnification Module Fields of View for Various Flip Lens Magnification Factors 


Smallest 
Middle ? 
Largest 


: £ismetg[ 
12 mm 
20 mm 
30 mm . 


; j .4x :: g 

15.7 mm 

22 mm 

30.8 mm 


14.7 mm 
22 mm 
33 am ! 


13.75 mm 
22 mm 
35.2 mm ! 


1.7 X 

12.9 mm 
22 mm 
37.4 mm 


Ti ffing of the different sized chambers into a fixed size pupil is accomplished by running the beam backwards 
through the system. The Keplerian telescope, without the flip lens, relays an image of the surface in the middle 
•ized chamber. Inserting the flip lens, the system relays images of the larger or smaller sized chambers. For 
proper imaging the surface must be in the front focal plane of the lens closest to the oil. The resultant image 
will be in the bacV focal plane of the last lens of the projection/imaging module. This image is transferred to 
the camera mod ule by the pupil relay and measurement modules - forming a Ronchigram. 



5.5 Pupil Relay Module 


Tht miDil relay module is a unit magnification AFOCAL system. Its primary function is to extend the beam 
!£5! ^thTmlurement module is at a location where the mission specials Rinsed t the gratrng. 

Additionally, it serves as an alternative measurement module J ho i e deformations, 

the beam not accessible with the measurement module. This is required to prope y 


5.6 Measurement Module 

chances the equivalent wavelength. The equivalent wavelength depends only upon the grating e spacrng 

grating murtbe pl^at a greater distance from the measurement module focus than a fine grating, 
distance from focus for a given grating depends entirely upon its tine spacing: 

I = 10 d {FI#) ■= 22 JSd 

Table 5 summari zes the distances needed for each grating to produce 10 lines. 

Table 5. Distance from Focus for Spatial Sample Points 



To reconstruct a two-dimensional surface contour, separate horizontal and verticd gratings, or crossed gratings 
can be inserted into the beam. Radial or circular zone gratings may be tned. 

5.6.1 Fill Level Indication 

— -.Ear isss °f of 

to wrtoe will b. to. If <!» «*« wU! f< T * ‘ S^aJisSXib « Mission 

Fa. of to — to totoo. ? '•£ ‘^oT^Z . to to 
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under or overfilled, the lines will reappear and increase in number. 

























5.7 Camera Module 


The Ronchigram produced by the measurement module is imaged onto a diffuser plate where a real image is 
formed. This image is viewed by a video camera and recorded on video tape for subsequent (post-flight) data 
analysis \ The physical distance from the diffuser screen to the camera is determined by the focal length of the 
ouoeni Jens, the camera sensor format, and the size of the image on the diffuser screen. 


d 


1 ♦ 



For the flight system, the diffusa* image radius is 10 mm, the tensor radius is 2.2 mm, and the camera lens 
focal length is 28 tm Thus, the stand-off distance is 155 mm. The camera is a 1/2 inch format RS-170 video 
camera. 

In addition to bang able to spatially resolve the surface feature without vignetting, the system is required to 
temporally resolve the oscillation. It is assumed that the period of this oscillation is on the order of 2 to 5 
seconds. Thus, a standard video camera with a 60 Hz field rate is adequate. Finally, it is assumed that there 
may be high-speed temporal gitter in addition to the lower-speed temporal oscillation. This gitter could reduce 
the contrast of the pattern. To eliminate gitter the camera will need a high speed shutter and the source will 
need sufficient output power to expose the camera. 

6.0 CALIBRATION 

If the optical system were aberration free, the Ronchigram of an ideal flat mirror would be a series of perfect 
straight lines. However, given the total number of positive optical components in the optical system, such a 
pattern is impossible. By measuring the beam reflected from a ’perfect’ reference-flat insert, after assembly, 
these e r ro r s can be characterized and removed from all measured data. Additionally, the system can be 
calibrated by testing a senes of known diamond-turned test plate samples. 

Calibration of the prototype instrument was accomplished by measuring a known amount of defocus aberration 
introduced by translating the collimating leas: 


_i — 

8 (F/#) 1 

Longitudinal Aberration 

The result of this calibration gives a measurement accuracy of approximately ±10%. Measurements made on 
the known test samples confirm this result. 

He system is designed for die insertion of interchangeable gratings at various locations in the beam. Thus, 
grating alignment errors can be a problem. To assess this, the data analysis process has been evaluated for its 
sensitivity to errors in the position of one grating relative to its orthogonal partner. In general, the Ronchi test 
is relatively insensitive to small grating placement errors (except for translation along the surface normal ) for 
grating frequencies less than 3 to 5 Ip/mm. 


W, 


020 


where: c x ■ 


7.0 CONCLUSIONS 

A prototype Ronchi instrument to monitor steady-state and oscillatory deformations, and fluid fill level of a free- 
surface of oil has been developed, demonstrated, and characterized for STDCE-2. A single-channel Ronchi 
configuration has been selected as the flight instrument. The flight instrument is currently being designed. 



n,e Roochi test is • shearing interferometer. The number of sheared wavefronts and the magnitude of the shear 
is dependant upon the characteristics of the grating. The fringe spacing (equivalent wavelength) u» ■ 
R 0 X«m Unproportional to the spring of the grating and the F/f of the test beam. The dynamic range and 
■ensitivfrTof die Roochi test can be changed by simply changing the grating spacing. The pro IT* 
SSS^chi system has an accuracy of ± 5%. For a shewing ^stem, the nnoumum memorable error is 

determined by the surface curvature. 
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Figure 1. Surface Tension Driven Convection Experiment. 
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Figure 2. Comparison of STDCE-1 and STDCE-2. 
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Figure 3. Single Channel Ronchi Test. 


Figure 4. Two Channel Ronchi Test. 
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Figure 5. The Ronchi Test 



Measurement of a 25 micrometer bump with three different gratings at a fixed distance z from 
focus producing three different equivalent wavelengths. 
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Figure 8. Sche mati c Diagram of the Illumination Module. 
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Figure 9. Schematic Diagram of the Projection/Imaging Module Magnification Optics. 
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Figure 10. Schematic Diagram of the Polarization Module Component Placement. 





Contouring of a free oil surface 
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ABSTRACT 

A non-contact (Le. optical) technique is required to measure the deformation of a free oil surface 
under the influence of a localized thermal load. This deformation is caused by surface tension driven thermal 
convective flow inside the fluid, and can be as large as 250 micrometers. Therefore, conventional interferometry 
is not possible. Instead, a Ronchi technique is proposed for contour mapping the oil surface. This paper presents 
a two channel Ronchi instrument, and some preliminary results from the analysis of two aspheric surfaces, a diamond 
machined metal surface and the free oil surface. 


L INTRODUCTION 

Scientists at Case Western Reserve University and NASA Lewis Research Center are conducting research 
on a variety of fundamental phenomena inherent in low gravity processes. One phenomenon being studied is 
thermocapillary flow, or Quid flow driven by surface tension. Thermocapillary flows result when temperature 
differences arise across a liquid/gas interface. These differences cause a change in the surface tension along the 
interface, which results in thermocapillary Dow. Improvements in the understanding of thermocapillary flow may 
lead to technological breakthroughs in the space processing of semiconductor crystals. These experiments were a 
part of the 13-day Microgravity Laboratory (USML-1) shuttle mission in the summer of 1992. The experiment 
described in this paper was not on that mission, but may be included on future shuttle missions. 

One of the parameters of thermocapillary flow to be studied is the ability of the free surface to deform 
during flow. The development and implementation of a non-contact (Le. optical) method to monitor these surface 
deformations as a function of experimental parameters such as fluid properties, heating modes, is crucial to the 
understanding of the physics of thermocapillary flow. 

The experimental surface is a three centimeter diameter pool of silicon oil. Deformations as large as 250 
micrometers are produced on the surface of the oil by thermocapillary flow. An optical technique instead of a 
mechanical technique was chosen to contour the free surface because touching the surface with a mechanical probe 
would deform the surface even more. There are two basic categories of optical techniques: geometric and 
interferometric. A geometric method measures surface slope, while an interferometric method measures surface 
height The geometric approach was chosen over the interferometric technique for two reasons. First, conventional 
visible wavelength interferometry would produce a fringe pattern too dense to analyze because of the anticipated 
deformation range of up to 250 micrometers. And, while infrared interferometric techniques could produce a 
resolvable fringe pattern, the infrared beam might affect the system’s thermal loading. Second, geometric techniques 
are dependant on the slope of the surface and are therefore insensitive to any jitter or tilt experienced by the oil 
surface while in space. One disadvantage of using geometric techniques is that two orthogonal data sets, either x- 
slope and y-slope or r-slope and ©-slope are needed to completely reconstruct the surface shape. 

Five c las** 1 ** 1 geometric techniques were considered: Foucault knife-edge test, Hartmann screen test, wire 
tests, grating interferometry, and Ronchi test The knife-edge and wire tests were not chosen because they provide 
data for only a single zone or diameter, not for the entire aperture. The Hartmann screen can provide two orthogonal 
data sets at one time, but it has a limited spatial resolution when compared with the Roochi test. The Rondutest 
was chosen as the technique to use for three reasons: the small dimensions of the experimental pool are tailored to 
the Ronchi requirement of a collimated beam, the Rchchi data is less complicated to analyze, and it can be integral 
with phase-measuring techniques. 
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2. BACKGROUND 


2*1 Relevant work 


The Ronchi test was discovered by Italian physic! 
Vasco Ronchi in 1921. Ronchi was attempting to measure 
Ihe radius of curvature of a concave spherical mirror. Upon 
looking through the grating at (he mirror, the surface defects 
of the mirror became apparent. The mirror surface appeared 
lined with twisted and irregular fringes. 1 Since its 
conception, the Ronchi test has been applied to die 
measurement of (he aberrations or defects in many different 
applications and using many different test setups. 

A single channel Ronchi test is shown in Fig. 1. 

A collimated wavefront is incident upon a test surface. The 

reflected wavefront contains information about die surface 
contour. The collecting lens focuses the wavefront onto a 
Ronchi paling. The irradiance pattern produced on the 
observation screen is a sheared interferogram called a 
Ronchigram, where the irradiance across the pattern is 
proportional to the slope of the test surface in the direction 
perpendicular to the grating structure. 


CotllMtlng Lane 




Figure 1 Single channel Ronchi test 


3. RONCHI THEORY 


3.1 Linear system theory 


The Ronchi test can be theoretically explained as either a geometric ray or an interferometric test. From 
ageometnc ray point of view, the fringes in a Ronchigram are the result of the deviation of a ray from its ideal path 
because of slope errors. In an interferometric test the fringes in the Ronchigram are produced from the interference 
betwtren overlapping diffraction orders. The overlapping wavefronts are produced by a shearing of the original 
wavefront Therefore the Ronchi test can be approached as a sheared wavefront interferometer. 



The basic function of the Ronchi test can be explained from a linear system theory point of view using the 
geometry shown m Fig. 2. The wavefront reflected off of the test surface. u(x,.y t ). contains information about the 
surface contour. The wavefront observed at the focus of the lens, directly before the grating, is the Fourier transform 
of Inc input wavefront. 




(i) 




After the grating, the wavefront is modulated by the grating transmission function, t(x2.y2). 
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Al the observation screen, the wavefront amplitude is the Fourier transform of the wavefront at the focal 
plane of the lens. 
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The wavefront at the observation screen is a scaled version of the incident wavefront, u, , convolved with 
the Fourier transform of the grating transmission function. 


«(j,^ = «< ^ ) . nj-Y-) 

z z Xz Xz 


(3) 


When the distance z ( distance from the grating to the observation screen ) is equal to the focal length of 
the collection lens, the wavefront scale factor disappears and equation 3 is the original wavefront convolved with 
the Fourier transform of the grating transmission function. The grating transmission function is the important 
mechanism in the Ronchi test 


The irradiance at the observation screen is given by: 

E&yy j) * < I«(Wj)I i> 


(4) 


3 2 Square wave grating 


The type of grating used in the Ronchi test determines 
A square wave grating was used in this experiment A square 
straight line bands. The Fourier transform of a square wave 
function. The square wave function and its Fourier transform 

When the width of the bands in the square 
wave grating is equal to ope half of the grating period, 
the even orders of the sine function are missing. If the 
contributions of the higher odd orders ( 3, 3, 7, ..) are 


the interpretation of the fringes in the Ronchigram. 
wave grating is a series of transparent and opaque 
grating is a sine function convolved with a comb 
ate shown in Fig. 3. 



Figure 3 50/50 Duty cycle square wave grating 



assumed to be negligible, the Fourier transform of a square wave grating can be approximated as three delta 
functions. Plugging this result into equation 3 gives the wavefront at the observation screen for a square wave 
grating with a 50/50 duty cycle. 


«*(*,. 7j) ■ *0^) * «(^ [Xj - <fl) * [x, ♦ d]) 

Equation 5 is of the same form as the general equation for a lateral shearing interferometer. 1 

"I* - y. y) - "T* * j. y) - mX 


( 5 ) 


( 6 ) 


The A in equation 6. referred to as the shear distance, is the distance between the centers of the wavefronts. 
Figure 4 shows a schematic of the shearing effect of the square wave grating. 


U3 



The iiradiance at the observation screen is given in equation 1. Just as in a shearing interferometer, the 
interference fringes in the Ronchigram are produced in the overlap region of the laterally sheared wavefronts. These 
overlapping regions are expressed in the irradiance equation as the cross terms in parenthesis. 

E(Xj, y 3 ) * E M + E k + E t + (uji 3 * u‘uj * {ujt’ + u’uj * («*n c * ♦ «*\) ^ 


There are three fringe patterns produced in the Ronchigram, the interference between wavefronts (a4>). (ax), and 
(b,c). Given our setups coarse gratings (50 lp/in), the actual shear distance is very smalL The observed irradiance 
(Ronchigram) is produced from the overlap region of all three wavefronts. 

3 3 Equivalent wavelength 

In interferometry, the fringe spacing in an interferogram is proportional to die wavelength of light used in 
the interferometer. The irradiance at any point in an interferogram, E(x.y). can be expressed as a function of 
measurement wavelength and optical path difference (OPD) between die two interfering wavefronts. 

EOcj) -EJ 1 ♦ cos( - - g— )) (*) 

where X * measurement wavelength 


The OPD in a shearing interferometer is given by the difference between the wavefront and the sheared wavefront. 

OPD - JUXy) - Wfr ♦ AjyO - AWXy) W 

where A * shear distance 

The relationship between the differential wavefront and the slope of the wavefront is shown in equation 10. 


AWXrj) - ISiS& Ax « «(xj) Ax (10) 

ox 

where a(xj) • wavefront slope 


Inserting equation 10 into equation 8 gives the inadiance equation for a shear interferometer. This is also 
the inadiance equation of a Ronchi test, since the Rone hi test is similar to a shearing interferometer. 




Ejtl ♦ cos(^f^)) 


( 11 ) 


The equivalent wavelength, of a shearing interferometer is defined as the measurement wavelength 
divided by the shear distance. 


. X (12) 

X *= — 
m Ax 


Inserting the equivalent wavelength definition into equation 1 1, the inadiance equation can be rewritten as equation 
13. 




E,( 1 * cos(^)) 


(13) 


Thus, by proper scaling of the interfeiogram, conventional interferogram analysis techniques will provide a direct 
measurement of the wavefront slope. 


33.1 Equivalent wavelength (Square wave grating) 

In order to derive an equivalent wavelength for a square wave grating with a 50/50 duty cycle the dominant 
fringe pattern in the Ronchi gram has to be determined. The shear between wavefronts (a,b) and (a,c), in Fig. 4, is 
dominant over the shear between wavefronts (b,c) for two reasons. First, the intensity of the wavefront u,, the zero 
diffraction order, is much higher than the intensity of the other two wavefronts. Second, the fringe pattern produced 
from (aj>) and (a,c) are identical and they overpower the fringe pattern produced from (b,c). 


The shear between wavefronts (a,b) and (a£) is identical and is equal to 



The shear distance is normalized by dividing ihc shear by the beam radius. Normally, the beam radius is 
equal to the radius of the collection lens because the lens pupil is over filled and the distance z is equal to the focal 
length of the collection optics. 


a' Ajc 
Ax * 
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Substituting the normalized shear into equation 12, the equivalent wavelength for the Ronchigram produced from 
(a,b) and (a,c) is 


X - -J-. < 15) 

* 2m 

The fringe spacing in the Ronchigram is independent of the measurement wavelength. It is only dependant 
upon the spacing of the grating and the f/# of the test beam. This means that the Ronchi test can be performed with 
either a monochromatic or white light source. In addition, the sensitivity of the Ronchi test can be changed by 
simply changing the grating frequency. 


4. TWO CHANNEL DESIGN 

In order to reconstruct the surface contour, two orthogonal slope patterns, either x -slope and y-slope or r- 
slope and 8-slope, are required. A single channel Ronchi test can be used to characterize a static surface by rotating 
the grating 900 during the experiment to obtain two orthogonal slope patterns. The oil surface is a dynamic surface 
and requires that two orthogonal slope patterns be obtained in real time. A two channel Ronchi test has two separate 
arms, one with a vertical grating and one with a horizontal grating. The two channel lest can produce two orthogonal 
slope patterns in real time. 

4.1 Optical layout 

A schematic of the two channel Ronchi test is shown in Fig. 5. A collimated wavefront is incident on the 
oil surface. The wavefront reflected off of the surface contains information about the surface contour. The reflected 
wavefront is split and directed through two separate channels. One of the wavefronts is focused onto a horizontal 
grating by collecting lens, L3, producing a Ronchigram with vertical slope information. The other wavefront is 
focused onto a vertical grating by collecting lens, L2, producing a Ronchigram with horizontal slope information. 
The two channels are staggered so that the distance from the oil surface to each collection lens is equal. The 
wavefronts leaving lenses L4 and L5 are positioned side-by-side on an observation screen. 



Figure 5 Two channel Ronchi test 





4 2 System tltt 


The surface under test should be located in the front focal plane of the collection optics. The maximum 
surface slope that will enter the collection lens of the Ronchi test is dependant upon the W of the collection lens. 
From the geometry shown in Fig. 6, it is clear that the f/H requirement to sec a given slope is. 

L . 1 ( 1<0 

t 2 tan(2a) 

A graph of the system 17# required for a given surface slope is provided in Fig. 7. The two channel Ronchi test 
designed for this experiment is a f/4 system. Therefore a surface slope of up to approximately 60 pm/mm will enter 
the two channel Ronchi test The maximum expected slope of the oil surface is 26.7 pmAnm. 



Figure 6 System geometry 


F-number Requirements 

"for surface slope" 



43 Grating frequency selection 

One of the advantages of the Ronchi test is that the dynamic range and resolution of the system are 
dependant on the characteristics of the gratings and the 1/# of the system not the measurement wavelength. Therefore 
a simple change of the gratings will permit a different range of deformations to be monitored. This gives the Ronchi 
system great flexibility. 

The proper grating spacing, d, for a given maximum surface slope can be found using equation 17. To 
enhance data evaluation, the optimum number of fringes per pupil radius is 10 fringes. Fewer fringes would reduce 
the accuracy and more fringes may introduce additional complexity. 

a- 2/g ~ 
m 

where f "focal length, a^, • maximum surface slope and m * fringes per r adius 


( 17 ) 




5. TWO CHANNEL CALIBRATION 


5.1 Calibration experiment 

The two channel Ronchi test was calibrated by testing a known amount of introduced defocus. The defocus 
was introduced by translating the collimating lens, LI, in Fig. 5. Defocus was chosen to calibrate the system for 
three reasons: 1) It is easy to introduce. 2) It is easy to analyze. 3) The amount of defocus introduced is directly 
proportional to the magnitude of the translation. Three different frequency square wave grating were used in the 
experiment Each of the different gratings has a different equivalent wavelength and will therefore provide a test 
of the equivalent wavelength equation for square wave gratings with a 50/50 duty cycle. Table 1 is a listing of all 
the trials ran in the calibration experiment 


1 

Table 1 Listing of Trials la Two Quoad Calibration Experiment 


Trial 

Grating Frequency 

Equivalent Wavelength 

Translation 

Deiocos 


ft/®] 

[micrometers] 

M 

[Waves] 

I 

250 

12.10 

2-0 

1.045 

2 

250 

12.10 

3.0 

1.623 

3 

250 

12.10 

4.0 

1221 

4 

250 

12.10 

5.0 

2.903 

5 

150 

20.16 

3-0 

0.950 

6 

150 

20.16 

4.5 

1.509 

7 

150 

20.16 

&0 

2.049 

8 

150 

20.16 

7 JS 

2^88 

9 

50 

60.48 

5-5 

0.650 

10 

50 

60.48 

10 

0.882 

n 

50 

60.48 

9-5 

1.076 


52 Ronchigram analysis 


The two channel Ronchi test produced two Ronchigrams far each trial. These Ronchigrams were of 
orthogonal slope patterns, x-slope and y-slope. Each slope map was digitized and fit to a set of Zemike polynomials. 
The two slope maps were then combined to produce the original surface contour by linear combinations of the 
Zemike polynomials. 3 

5J Calibration results 

Figure 8 is a graph of the data shown in Table 1. The trials for each grating frequency demonstrate a linear 
relationship between the amount of lens translation and the amount of introduced defocus. Multiplying each defocus 
value by the appropriate equivalent wavelength value produces the result shown in Fig. 9. All the trials fall upon 
the same line. The equivalent wavelength equation does relate data taken with different grating frequencies. 

Figure 10 shows the percent error between the data and the best fit line shown in Fig. 9. The average 
percent error of the current two channel system is ±5%. The percent error for small amounts of defocus is higher 
than five percent. It is speculated that this is because of the inaccuracy involved in placing the gratings at best focus. 






























































Two Channel Ronchi Calibration Two Channel Ronchi Calibration 


Defocus Scaled by Equivalent Wavelength 
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Figure 8 Two channel calibration data Figure 9 Data scaled oy equivalent wavelength 


Two Channel Ronchi Calibration 

Percent Error 



Figure 10 Percent error 

6. PRELIMINARY DATA 

The equivalent wavelength equation needs to be modified when the test object is a surface instead of a 
wavefront as in the case of the two channel calibration. The cosine of the reflected angle is multiplied by the 
denominator of the equivalent wavelength equation. When the reflected beam is normal to the test surface equation 
IS is the proper scale factor (cos (0) * 1). 


X - 20Wcos(e) 

where 0 * the angle between the optic axis of the collection lens and the surface normal 


( 18 ) 





6.1 Diamond machined known -surface 


A diamond machined aluminum surface was examined with (he (wo channel Ronchi (est. The surface was 
machined flat except for a 5 mm diameter area in the center of the surface. This central area contained a 54 
micrometer depression. A Ronchigram of the surface and a contour map of the surface generated from the 
Ronchigram information are shown in fig. 11. The two channel system measured the depression to be 52.412.6 
micrometers. The area around the depression is flak 


6 2 Ofl surface 


The oil surface usually has either a depression or a bump at the location the COj beam strikes the surface. 
The size of the feature can be changed as well as the polarity by adjusting the size of the incident CO, beam and 
the laser power level, fig. 12 shows the Ronchigrams and a contour map generated from the Ronchi grams for a 
bump on the surface. The COj beam size on the surface was 6 mm. The two channel system measured the bump 
to be 1711 micrometers. 




Figure 11 Diamond turned surface 


Figure 12 Free oil surface 




7. CONCLUSIONS 


The Ronchi test can be treated as a shearing interferometer. The number of sheared wavefronts and the 
magnitude of the shear is dependant upon the characteristics of the grating. The fringe spacing in a Ronchigram 
is proportional to the equivalent wavelength of the Ronchi system. The equivalent wavelength is dependant upon 
the spacing of the grating and the f/# of the test beam. The sensitivity of the Ronchi test is both known and variable. 
It can be changed by simply changing the grating spacing. 

The two channel Ronchi system was calibrated and has an accuracy of ±5%. The preliminary data 
demonstrates the ability of the Ronchi lest to contour an aspheric surface. 
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ABSTRACT 

A noo-cootact 0-e. optical) technique is required to measure the deformation of a free oil surface under the 
influence of a localized thermal load. This deformation is caused by surface tension driven thermal convective flow 
i nri*. the fluid, and can be as large as 250 um. Therefore, conventional interferometry is not possible. Instead, a 
Ronchi technique is proposed for contour mapping the oil surface. This paper presents a design coocept for a two 
channel Ronchi instrument and some preliminary results from a feasibility study using a single channel Roochi 
instrument 


1. INTRODUCTION 

Scientists and engineers at Case Western Reserve University and NASA Lewis Research Center are 
conducting research on a variety of fundamental phenomena inherent in low gravity processes. An important 
phenomenon considered to have a significant impact on processes such as containeriess processing and crystal growth 
is that of thennocapillary flow, or fluid flow driven by surface tensioo. Thennocapillary flows result when 
temperature differences arise across a liquid/gas interface. These differences cause a change in the surface tension 
along the interface, which results in thermocapillary flow. Improvements in the understanding of tbennocapillary 
flow will help lead to technological breakthroughs in space processing of semiconductor materials. 

The ability of the free surface to deform during thennocapillary flow is an important aspect of this problem. 
These deformations can alter the surface temperature distribution, the driving force in the flow, and the resulting fluid 
flow patterns , causing an unsteady flow which is detrimental to the crystal growth process. Therefore, the 
development and implementation of a technique to study the surface deformation as a function of experimental 
parameters such as fluid properties, beating modes, is crucial to the understanding of the underlying physics. 

This paper describes the design of a technique for measuring surface deformations of a 3 cm diameter pool 
of silicon oil. Since the surface is a fluid, touching it would change its shape. Therefore the technique must be non- 
contact, i.e. optical. Two techniques were considered: geometric and interferometric. A geometrical technique 
measures surface slope, while an interferometric technique measures surface height The geometric approach was 
chosen over the interferometric approach for two reasons. First, conventional visible wavelength interferometry 
would produce a fringe pattern too dense to analyze because of the anticipated deformation range of up to 250 
micrometers. And, while infrared interferometric techniques could produce a resolvable fringe pattern, the infrared 
beam might affect the system’s thermal loading. Second, geometric techniques are dependant on the slope of the 
surface and are therefore insensitive to any jitter or tilt experienced by the oil surface while in space. One 
disadvantage of using geometric techniques is that two orthogonal data sets, either x -slope and y-slope or r-slope and 
6-slope are needed to completely reconstruct the surface shape. 

Five classical geometric techniques were considered: Foucault knife-edge test, Hartmann screen test, wire 
teste, grating interferometry, and Ronchi test The knife-edge and wire teste were not chosen because they provide 
data for only a single zone or diameter, not for the entire aperture. The Hartmann screen can provide two orthogonal 
data sets at one time, but it has a limited spatial resolution when compared with the Roochi test The Ronchi test 
was chosen as the technique to use for three reasons: the small dimensions of die experimental pool are tailored to 
the Roochi requirement of a collimated beam, the Roochi data is less complicated to analyze, and it can be integrated 
with phase-measuring techniques. 
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2. BACKGROUND 


The Ronchi test was discovered by the Italian physicist Vasco Roochi in 1923. He noticed that when a 
grating was positioned near the center of curvature of a minor that the image of the grating was superimposed upon 
the original grating. The combination of the two gratings produced a moire type effect whose pattern depended upon 
the aberrations of the mirror. Roochi felt that this method was a good check for mirror quality. 1 

The properties of the Roochi test can be 
understood by examining the single channel Roochi 
setup shown in Fig. I. The source in Pig. 1 is a 'point 
source' created by spatially filtering a HeNe (X — 632.8 
nm) laser beam. A lens is positioned ooe focal length 
in front of the spatial filter, producing a collimated 
h>nm of light The collimated wavefront is incident 
upon the test surface and reflected to an achromatic 
imagi ng lens. The reflected wavefront contains 
information about the surface deformations. The 
imagin g lens focuses this wavefront onto a Ronchi 
grating. The irradiance pattern produced on the 
observation screen is a sheared interferogram called a 
Roochi gram, where the irradiance across the pattern is 
proportional to the slope across the test surface. 



The Ronchi test measures wavefront aberrations present in a system. Therefore, if there are no aberrations, 
the Ronchigram will have no fringes. When the lest surface in Fig. 1 is a plane minor no aberrations are introduced 
into the collimated wavefront Therefore, if the grating is at the focus of the imaging lens, either a bright spot or 
a dark spot is produced on the observation screen. A dark spot indica t es that the focused beam is incident upon a 
dark ba nd in the grating and a bright spot indicates that the focused beam is incident upon an area between two dark 
fa»nHc in the grating. Now, if the grating is moved away from the focal plane of the imaging lens, defocus is 
introduced. Defocus appears as straight line fringes in a Ronchigram. 

When the surface is aberrated, the Ronchigram contains fringes which are proportional to the slope of the 
surface. Since the Roochi test measures the wavefront’s slope in a direction perpendicular to the grating lines, it 
is necessary to obtain two orthogonal patterns to fully represent the surface shape. 

3, CERTIFICATION OF A SINGLE CHANNEL RONCHI TEST 

For any measuring instrument, it is necessary to develop a certification technique. For the Ronchi test, this 
can be accomplished by introducing a known amount of aberration and then testing to see if that amount is actually 
measured. In this section, a certification process using defocus is discussed. 


3.1 Theory 


An ideal spherical wavefront is ooe which will produce a perfect point image at its focal point Such a 
wavefront is described by, 


w(x,y) - 


x 2 * y 2 
2R 


where (x,y) are pupil coordinates, and 
R is the wavefront radius of curvature. 


( 1 ) 




If the wavefront is defocused by an amount Aw(x,y), 

( 2 ) 

w'(x,y) - w(x,y) * Aw(x,y) 

it will focus at a slightly different location, R+e,, where c, is called the longitudinal ray aberration. The wavefront 
error introduced by defocus can be described by, 

Aw(.xj) ■ (i 1 ♦ y*) ^ 

where w^ is the Defocus coefficient. 

Figure 2 shows a reference wavefront and an aberrated wavefront with a profile function of W(x,y). The longitudinal 
aberration, e,, and the transverse aberration, e^, are also shown. The slope error, a, is the angle between the 
intersection of the reference wavefront and the aberrated wavefront surface normals. The slope error is equal to the 
partial derivative of the profile function W with respect to y. 

a w <0 


W(x,y) 


Reference Wavefront 
-Aberrated Wavefront 


y 

\zl 


r - pupil radius R - radius of curvature 
\ x, y - pupil coordinates 

,,x ez - longitudinal aberration 

ey - transverse aberration 
a - slope error 

Figure 2 Comparison of reference and aberrated wavefronts 

Using similar triangles wd die small angle approximation it can be shown from Fig. 2 that: 

ft 

e. - o R e. - — e 
y 1 r 7 

where r * pupil radius. 


( 6 ) 


Since Roochi grams are in normalized pupil coordinates, it is necessary to normalize x, y, and cl 


, x , y 1 bW 

x' ~ — y' - L a - - — — 
r r r dy 


where r « pupil radius and (x'y 1 ) - normalized pupil coordinates. 

By substitution of equations 4 and 6 into 5, the longitudinal and transverse aberrations can be rewritten as: 


„ dW 

R 2 W 
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(7) 


( 8 ) 


Taking the partial derivative with respect to y of the defocus term in equation 3, writing it in normalized pupil 
coordinates, and substituting it into equation 8 yields equation 9. 

(9) 


pi 

'.-2^7 


The f-number of the system, (f/#), shown in Fig. 2 is the radius of curvature of the reference sphere divided by the 
pupil diameter (2r). 

/ . ± m 

♦ ' 2r 


Substituting the f-number into equation 9 provides an equation for the defocus in a wavefront that is dependant only 
upon the amo unt of longitudinal aberration and the system f-number. 
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Thus, by dimply translating the Roochi grating, a known amount of aberration can be introduced. 


33. Measuring defocus 


The number of fringes in a Roochi gram is related to the partial derivative of the wavefront aberration 5 . 


bW m _ md 

dy R 

where R - radius of curvature m - number of fringes per radius, 
and d ■ grating spacing. 


( 12 ) 


Writing equation 12 normalized pupil coordinates gives: 

dW 


m d r 


(13) 



Taking the partial derivative of the defocus term in equation 3 and writing it in normalized coordinates yields: 


3Aw(x,y) 

3y' 


2 w. 


(00 


( 14 ) 


Setting the right rides of equations 13 and 14 equal and substituting in the definition of the W (equation 10) results 
in the following relation: 




m d 

4 (L) 

# 


(15) 


Thus, the amo unt of defocus in the wavefront can be dete rmined by counting the number of fringes per pupil radius 
in the observed Roochigram. 

33 Certification experiment 


The angle channel Roochi test shown in Fig. 
1 was certified by testing a known amount of 
introduced defocus. The grating was initially placed in 
the focal plane of the lens and its position was recorded 
for a reference value. The grating was then displaced 
a known ■mnnnt towards the lens and the output fringe 
pattern was photographed. An example defocus fringe 
pattern is shown in Fig. 3. Tran sl ating the grating is 
equivalent to introducing longitudinal aberration, 
into the system. The number of fringes can be 
tetemuaed directly from the photographs for each 
t ransl ation The experiment was performed with 
different frequency gratings and different grating 
orientations to insure' that equation 15 was correct and 
to get a visual feel for the effects of using different 
gratings. 

Table 1 presents the experimental data for 
various translations and grating frequencies. Equation 
Figure 3 Defocus fringe pattern jj was use( j to calculate the value of defocus from a 

known amount of translation. Equation 15 was used to 
calculate the value of defocus from an observed number of fringes. The difference between the two calculations is 
shown in the defocus difference column in Table I. The reason that there are two different values for the f-number 
is because a portion of the data was taken with the grating in the vertical direction and a portion was taken with the 
grating in the horizontal direction. The elliptical beam produced by the test surface, which is oriented at 45°, results 
in two different f-oumbers’s for the different grating orientations. 

The defocus values from the two different methods compare very well. The largest defocus difference is 
38 microns. The orientation of the grating and the grating frequency do not aeem to have any effect on the accuracy 
of calculated defocus values. The results show that the theory is correct and that the amount of defocus can be 
accurately calculated 




Table 1 Comparison of Calculated Defocus Values 


Trial 

Orating 

Spacing 

(Unes/mm) 

F-numbcr 

2.17*0.05 

3.08±0.09 

Translation 

(*0.o1* mm) 

Observed 

Fringes 

Defocus 

from 

Translation 

(mm) 

Defocus 

from 

Fringes 

Defocus 

Difference 

(mm) 

I 

9.843 

2.17 

0.71 

4.5 

huh 

■ 

0007 1 

2 

9.843 

2.17 

1.62 

u 

00*3 

oo<* 

0.001 J 

3 

9.843 

2.17 

3.60 

13.0 


sm 

HffK 

4 

9.843 

2,17 

4.92 

180 


0105 

0025 

5 

1.969 

2.17 

5.16 

4.0 

0137 

0.117 


6 

1.969 

2.17 

8.40 

4.5 

0.223 

0-19Q 


7 

1.969 

2-17 

7,42 

5.5 

0.197 


HJjjjjHi 

8 

1.969 

2.17 

1 

8,5 

sm 

H5H 

1, 

9 

5.906 

3.08 

1.14 

4.0 

mi 

0.027 

0.012 

10 

5.906 

3.08 

3.78 

7.0 

mmm 

I 

© 

0.002 

11 

5.906 

3.08 

6.23 

12.0 

0002 

0.082 

0.000 

mmmw U acaM= _ 

5.906 


7.83 

15.0 

0.103 

0.103 

0.000 


4. TWO CHANNEL DESIGN 


In order to fully characterize the surface, it is necessary to measure the surface slope in two orthogonal 
directions. The single channel Ronchi system can be used to obtain two orthogonal slope patterns by simply rotating 
the grating 90° during the experiment However this method only works for static surfaces like the plane mirror 
presently being used to calibrate the system. Since the pool of oil is a dynamic surface two orthogonal fringe 
patterns must be obtained simultaneously in real time. A two channel Roochi test can produce two orthogonal fringe 
patterns in real time. 



Figure 4 Two channel Ronchi test 


































A schemati c of a two channel Roochi test is sbowD in Fig. 4. The two channel test was designed to fit on 
a standard 2’ by 3’ optical breadboard. The space requirements of the optical component mounts are not shown in 
Fig. 4, but they were taken into account when parameters such as focal length of the lenses and the diameters of the 
elements were chosen. 

The source shown in Fig. 4 is a HeNe (k — 632.8 on) laser. The light emitted by the laser strikes a fold 
mirror (Ml) and is directed into a spatial filter (SP). The collimating lens, LI, is positioned one focal length away 
from the spatial filter, thereby producing a collimated beam of light The collimated beam or plane wavefront is 
incident upon the test sample at 45°. The wavefront reflected off of the test sample contains information about the 
contour of the test surface. The reflected wavefroot is then divided into two wavefronts by the beamsplitter (BS) 
and directed through two separate arms. One of the wavefronts is focused onto a vertical grating by L2 producing 
a Roochigram with horizontal slope information, while the other wavefront is focused onto a horizontal grating by 
L3 producing a Roochigram with vertical slope information. The two dope maps are imaged by lenses L4 and L5 
onto the observation screen. The lenses L1-L5 were chosen to be acbromats to insure that minima] amounts of 
spherical aberration are introduced. 


The wavefronts exiting L4 and L5 are positioned side-by-side on an observation screen as pictured in Fig. 
4 by the fold mirrors, M3 and M4. Additional optics are required in order to form an image of die two Ronchi grams 
onto a video camera. 

4.2 Sample and sample mount 

The test sample to be used in the two channel Roochi setup was designed and manufactured by 
undergraduate mechanical engineers and applied optics students at Rose-Hulman. It was designed to simulate the 
free oil surface. Three criteria were used in the design of the test sample: the surface must be 3 cm in diameter, 
it must be capable of creating both odd and even deformations, and it must be large enough to insure dial stress 
fringes from its mounting do not overlap with experimental deformation fringes. The test surface is a circular piece 
of plexiglas mounted between two metal rings. The back metal ring is connected to a tip-tilt mount A black circle 
of 3 cm diame ter is centered on the plexiglas surface to fiducial ize die lest surface. 

Deformations in the test surface are created by 
adjusting the position of one of three threaded screws 
which are glued into the back side of the plexiglas. 
The three screws are held in place with a threaded nut 
which is attached to a metal plaie across the back of 
the camp le mount The metal plate is thick enough to 
insure that the plexiglas surface will deform before the 
metal plate. The middle screw is located at the center 
of die 3 cm test area and the other two screws are 
vertically displaced 1 cm on each side of die center 
screw. The screws can be translated forward or 
backwards to create an even deformation in the surface 
(a b»mp or a dimple). An odd deformation (a S- 
shaped deformation) can be created by a combination 
of screw adjustments, screwing the top screw back and 
the bottom screw forward or vice versa. The 
arrangement of the screws in the back of the plexiglas 
is shown in Fig. 5. 


Cross Section View 

Mount 


| Three Deformation 
Screws 


Metal Plate 


Plexiglas 
Test Surface 


Mount 


Figure 5 Arrangement of deformation screws 


4.3 Vignetting 

The next step is to bring the two designs, optical and mechanical, together to insure that no unwanted 
vignetting is introduced into the system. Vignetting occurs when the outer portion of the beam is partially blocked 
by the components or component mounts in the system. When designing an optical system, vignetting calculations 




are important to insure that the sizes of the optical compooeots meet the beam size requirements of the experiment 
The actual vignetting calculations for the two channel Roochi test are not provided, but some of the guidelines 
followed are provided. 

First when ordering optical components it is a good rule of thumb to allow 209b of the component diameter 
to be used for the mounting of the component Second, be aware of any surfaces in the design that are tilted relative 
to the beam path, like the test surface and the beamsplitter in the two channel Ronchi design. The beam that strikes 
the test surface needs to be large enough to overfill what appears to be an elliptical test surface due to the tilt of the 
surface and the clear aperture of the beamsplitter, which is reduced because of tilt, must be large enough to allow 
the reflected beam to pass through the mounted beamsplitter. Thud, be sure that none of the components on the 
table are blocking a beam that is propagating towards another component For example, the mounts for minors M3 
and M4 in Fig. 4 were designed to insure that the mount for M4 did not block the beam from M3 on its way to the 
observation screen. Fourth, be sure that the output beams can be positioned on the output device. In Fig. 4 the 
output device is the observation screen. 


4.4 Grating frequency selection 


One of the major advantages to using the Roochi test is that the dynamic range and resolution of the system 
are dependant only on characteristics of the grating used. Therefore if the experiment ever requires that a different 
range of deformations be monitored then only the grating needs to be changed, not the whole setup. This gives the 
system great flexibility. 


In order to choose the proper grating spacing ooe needs to know the maximum anticipated surface slope, 
the relative sphere radius of curvature, and the maximum number of fringes per pupil radius. T be proper grating 
spacing for a given experiment can be calculated using equation 16, where R is the reference sphere radius of 
curvature, is the maximum anticipated surface slope, and m is the desired number fringes per pupil radius. 


w 2Ra ~ 

a *= — 

m 


( 16 ) 


In the Ronchi test, the reference sphere 
radius of curvature is determined by the imag ing 
lens focal length. For the free oil surface, the 
anticipated maximum slope was found from Fig. 6 
to be 26.7 nmAnm To enhance data evaluation, 
the optimum number of fringes per pupil radius is 
10. Fewer fringes would reduce accuracy and 
more fringes may introduce additional complexity. 

4 £ Diffraction effects 


The Ronchi test can be viewed from a 
physical optics point of view as a diffraction 
phenomenon. Just as a diffraction pattern has 
higher diffracted orders so does a Roochigram. 
The clarity of the fringe pattern in the output 
Roochigram can be increased by removing some of 
the noise of the higher diffracted orders from the 
pattern. A spatial filtering experiment was 
performed to remove some of the noise and it had 
some positive results. A better approach is to use a sine grating instead of a step grating. The sine grating produces 
only the diffraction orders 0, and ±1. It will remove higher order noise from the fringe patterns and increase fringe 
clarity. 



Figure 6 Anticipated deformations 




5. CONCLUSIONS 


A two channel Roochi instrument has been designed to contour the free oil surface. It can produce two 
orthogonal fringe patterns at the same time, thus allowing a dynamic surface to be mooitored in real time. It is also 
designed to allow for easy implementation of a phase measuring capability in order to distinguish between bumps 
and dimples on the surface. 
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Sensitivity of two-channel Ronchi test to grating misalignment 
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ABSTRACT 

The Ronchi test is performed by placing a periodic grating at or near the focus of an optical 
system. Like most geometric tests, the Ronchi test measures wavefront slope instead of the 
contour of the wavefront like conventional interferometry. Therefore when a Ronchi test is used to 
determine the contour of a non-symmetric surface, two orthogonal slope interferograms or 
Ronchigrams are required to reconstruct the surface. A two-channel Ronchi test with a horizontal 
and a vertical grating can be used to reconstruct a non-symmetric test surface. This paper 
investigates the errors introduced by different rotational and translation misalignments of the 
gratings in the two-channel Ronchi test relative to one another. The Ronchigrams are analyzed 
with a modal technique based on Zernike polynomials. The severity of each misalignment is 
dependant upon the frequency of the gratings used and on the f/# of the system. 

I. INTRODUCTION 

The Ronchi test is performed by placing a periodic structure at or near the focus of an 
optical system and monitoring the interference pattern or Ronchigram produced. The Ronchi test 
can be explained as a geometric test or an interferometric test. From a geometric ray point of 
view, the fringes in a Ronchigram are the result of the deviation of a ray from its ideal path 
because of slope errors. As an interferometric test, the fringes in the Ronchigram are produced 
from the interference between overlapping diffraction orders. The overlapping diffraction orders 
are sheared versions of the original wavefront. Therefore the Ronchi test can be approached as a 
shearing interferometer. 

Shearing interferometers measure the slope of a wavefront instead of the height of the 
wavefront like conventional interferometers. One of the disadvantages of shearing interferometers 
is the requirement of two orthogonal slope maps of the surface, either x-slope and y-slope or r-slope 
and 6-slope, to completely reconstruct a surface without symmetry. If the surface under test is 
rotationally symmetric or if it has a static contour a single channel Ronchi test, Fig. 1, could be 
used to reconstruct the surface. The two orthogonal slope maps could be generated by rotating the 
grating ninety degrees during the experiment or by taking one photograph and analyzing it for 
both slope maps. However, if the surface under test is not rotationally symmetric and is dynamic, 
a two-channel Ronchi instrument must be implemented, Fig. 2. 

The main differences between the two-channel Ronchi instrument and the single channel 
instrument are the addition of a second measurement module and a pupil relay module. The relay 
module creates an intermediate image of the surface under test at the plane of the second 
beamsplitter. The image created by the relay module is the object for both of the measurement 
modules. The addition of the second measurement module allows orthogonal slope information to 
be obtained in real time. One of the measurement modules obtains a vertical, y-direction, slope 
map of the test surface, S, with a horizontal grating. The second module obtains a horizontal, x- 
direction, slope map with a vertical grating. The output of each channel is imaged onto an 
observation screen. It should be noted that the orientation of the image of the surface at the 
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Figure 1 Single channel Ronchi system 


Figure 2 Two-channel Ronchi system 


observation screen for the two-channel system is opposite to the orientation in the single channel 
due to the addition of the pupil relay module. 

The wavefront at the observation screen is the wavefront reflected off of the test surface 
convolved with the Fourier transform of the grating transmission function 1 : 

•Ow.) ■«<**) •*(&£) <U 

provided that the surface under test is in the front focal plane of the relay module and that the 
observation screen is in the back focal plane of the measurement modules. Therefore the grating 
transmission function is the important mechanism in the Ronchi test. It dictates how the fringes 
in the Ronchigrams should be interpreted. 

2. INTERFEROGRAM ANALYSIS 

Whenever two coherent wavefronts are at the same point in space, their amplitudes 
superimpose to create an interference pattern. The interference pattern contains information 
about the differences between the two wavefronts or their optical path difference (OPD). The two 
dimensional irradiance pattern created by the interfering wavefronts is known as an 
interferogram. 

2.1 Interferogram theory 

In interferometry, the fringe spacing in an interferogram is proportional to the wavelength 
of light used in the interferometer. The irradiance, E(xy), at any point in an interferogram can be 
expressed as a function of the measurement wavelength and the optical path difference (OPD) 
between the two interfering wavefronts. 






- E, (l ♦ co^^]] 
where X - measurement wavelength 


( 2 ) 


The OPD in a shearing interferometer is given by the difference between the wavefront and the 
sheared wavefront. 

OPD - WQcy) - W(?+Lxj) - A WQcy) 
where A - shear distance 


The relationship between the differential wavefront and the slope of the wavefront is given by: 

AB'i&y) - - a -ff ^ Aj - b(x^)Ajc 

dx 

where a(xy) ■ wavefront slope. 


Substituting equation 4 into equation 2 gives the irradiance equation for a shearing interferometer 
which iB also the irradiance equation for the Ronchi test. 
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where ■ — ■ equivalent wavelength 
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Thus, by proper scaling of the shear interferogram via the equivalent wavelength, conventional 
interferogram analysis techniques can be used to provide a direct measurement of the wavefront 
slope. 

2.2 Ronchi test equivalent wavelength 

To properly interpret a Ronchigram it is necessaiy to derive an equivalent wavelength, X*,, 
for the Ronchi test. The type of grating used in the Ronchi test governs how the fringes in the 
resultant interferogram or Ronchigram should be interpreted. An expression can be derived for 
the equivalent wavelength of the Ronchi system based on the characteristics of the grating and the 
geometry of the setup. 

Since the lateral placement of the grating in the Ronchi test is arbitrary, the amplitude 
sinusoidal grating can be modelled els Em amplitude cosine grating whose transmission function 
and Fourier transform are: 


4 


1 +C08 (^H)J - 6(x) ♦ 6(x-d) * 6(x*d) 

where d - the grating spacing. 


( 6 ) 



The Fourier transform of the amplitude cosine grating is a set of three delta functions. 
Substitution of equation 6 into equation 1 yields the expression for the wavefront at the 
observation screen. 
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( 7 ) 


The wavefront at the observation screen is three sheared versions of the original wavefront. The 
irradiance at the observation plane is given by: 




E(x 3 , yj m E a ♦ E b + E c + (u a Ut * k/u*) + («X + «.*»«) + K«c’ + «»«*)• 


( 8 ) 


The interference fringes in a Ronchigram are produced by the overlapping wavefronts 
which are the cross-terms in parenthesis in equation 8. The fringes correspond to changes in 
wavefront slope in a direction perpendicular to the shear direction. A low spatial frequency 
grating produces a small shear distance which in turn produces almost complete overlap between 
the wavefronts. High spatial frequency gratings produce large shear distances which reduce the 
amount of overlap. Therefore high frequency gratings do not test the full aperture of the 
wavefront. 

There are three sets of cross-terms in equation 8. Therefore there are three sets of sheared 
wavefronts: (a,b), (a,c), and (b,c). In order to derive an equivalent wavelength for the sine grating, 
the dominant fringe pattern in the Ronchigram must be determined. The shear between 
wavefronts (a,b) and (a,c) is dominant for two reasons. First, the intensity of the wavefront u t , the 
zero diffraction order, is much larger than the intensities of the other two wavefronts. Second, the 
fringe pattern produced from (a,b) and (a,c) are identical and they overpower the fringe pattern 
produced from (b,c). The equivalent wavelength is dominated by the shear between (a,b) and (a,c). 


The shear between wavefronts (a,b) and (a,c) is identical and is equal to : 
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where / - focal length, and d - grating spacing. 


( 9 ) 


Because interferograms are analyzed in normalized coordinates, the shear distance should be 
normalized by the maximum radius, the radius of the exit pupil of the system, which is the 
largest radius wavefront which can be tested. 
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Replacing x^ by D/2 where D is the diameter of the exit pupil of the system, gives: 
Substituting the normalized shear into the definition for X*, results in: 
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The equivalent wavelength of the Ronchi test is independent of the measurement 
wavelength. It only depends upon the spacing of the grating and the f/# of the system. This 
means that the Ronchi test can be performed with either a monochromatic or a white light source. 
In addition, the sensitivity of the Ronchi test can be changed by simply changing the grating 
frequency. 

2*3 Surface reconstruction from slope information 

Wavefront phase estimation from slope data is divided into two categories, modal and 
zonal. Zonal techniques estimate the phase value over a local region while modal techniques 
determine the coefficients of a set of aperture polynomials. 2 Rimmer developed a method for zonal 
estimation. 3 An example of a modal technique is the use of the Zernike polynomials to produce a 
global fit of the data. 4 ' 5 In both cases, a least squares estimation is used to reconstruct the 
wavefront phase. Work has also been published on the estimation of slope data with the use of 
Fourier transform techniques. 6 " 7 This paper uses a modal algorithm to reconstruct the wavefront 
shape through a linear combination of Zernike coefficients which are individually fit to orthogonal 
x- and y-slope data sets. 

All of the trials presented in this paper were analyzed with a modified version of the 
FRINGE software originally from the University of Arizona. Table 1 is a listing of the first eight 
Zernike polynomials. Table 2 provides the relationships between the Zernike coefficients of the 
slope data and the Zernike coefficients of the wavefront shape. For example, the wavefront x- 
astigmatism coefficient, A4, is dependant upon the x-tilt from the x-slope map, Bl, and the y-tilt 
from the y-slope map, C2. Finally, the Seidel coefficients can be calculated from the first eight 
Zernike coefficients as shown in Table 3. 

Table 1 First eight Zernike polynomials 


Term 

Aberration 

M 

0 

Piston 

l 

1 

x-tilt 

X 

2 

y-tilt 

y 

3 

defocus 

2(x 5 +y 2 )-l 

4 

x- astigmatism 

*v 

5 

y-astigmatism 

2xy 

6 

x-coma 

8(x*+xy*)-2x 

7 

y-coma 

8(x 2 y+y s )-2y 

8 

spherical 

6(x s +y s ) s -6(x a +y s )+ 1 



Table 2 Wavefront Zernike coefficients from linear combinations of slope Zernike coefficients 


Term 

Aberration 

Combination 

A1 

x-tilt 

cannot be determined 

A2 

y-tilt 

cannot be determined 

A3 

defocus 

(Bl+C2)/8-A8 

A4 

x-astigmatism 

(Bl-C2)/4 

A5 

y-astigmatism 

(B2-Cl)/4 

A6 

x-coma 

(B4+B3+C5)/9 

A7 

y-coma 

(B5-C4+C3)/9 

A8 

spherical 

(B6+C7)/16 


B coefficients are the x-slope. 

C coefficients are the y-slope. 

A coefficients are the surface shape. 

Table 3 Relationship between Seidel and Zernike coefficients 

Seidel term Zernike combination 

Tilt Magnitude ((A1-2A6) 2 +(A2-2A7) 2 ) W 

Angle Tan 1 ((A1-2A6)/(A2-2A7)) 

Defocus 2A3-6A8±((A4) 2 +(A5 ) 2 ) w 

The sign of the radical is chosen to minimize the magnitude. 

Astigmatism Magnitude ±2((A4) 2 +(A5) 2 ) iy2 

Angle 0.5 Tan 1 (A5/A4) 

The sign of the radical is opposite the sign used for defocus. 

Coma Magnitude 3((A6) 2 +(A7) 2 ) m 

Angle Tan* 1 (A7/A6) 

Spherical 6A8 

3. EFFECT OF GRATING MISALIGNMENT 

When aligned correctly the gratings in the two-channel Ronchi test should be at best focus of 
the system, perpendicular to the optic axis, and orthogonal to each other. However, misalignments 
can occur in the placement of the gratings which can produce measurement errors. The gratings 
have six degrees of freedom, three rotational and three translational as shown in Fig. 3. 

Experiments were performed with an f/4, two-channel system shown in Fig. 4 to 
investigate the effect of each misalignment on the system performance. In each case, the x-slope 
grating was used as the reference and the y-slope grating was manipulated relative to the 
reference grating. The gratings used in the experiments were 50 Ip/in, 50/50 duty cycle square 
wave gratings. The equivalent wavelength for these gratings can be approximated as the 
equivalent wavelength for the amplitude sinusoidal gratings. The collimating lens, Ll, was 
defocused for all the experiments so that there was a sufficient number of sample points across the 
aperture. This is analogous to displacing the grating away from focus. Defoe using the collimating 
lens produces an initial "tilt* fringe pattern to be produced in both channels. The number of tilt 
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Figure 3 Geometry to determine tilt fringes 



Figure 4 Two channel Ronchi system 


fringes can be determined from similar triangles in Fig, 5 and from a knowledge of the grating 
frequency used in the experiment. 

WjM - B1 - SL and ~C2-^ 

fix* J7# " rt4W ' ft# (12) 

where e x - the distance from the grating to focus, 
v » grating frequency, and fl# - the system /-number. 

Ub* 

^ The term is the initial aberration measured 

iy 2 in the x-slope or reference channel while W mHab)e 

| nxw is the initial aberration measured in the y-slope 

-Jc — or test channel of the two-channel Ronchi 

system. 

^ — o ^ 3d Rotation about the optic axis 

l j The first degree of freedom tested is a 

^ ^ rotation about the optic axis. The fixed x-slope 

Figure 5 Geometry to determine tilt fringes grating and the variable y-slope grating are at 

the focal plane of the system and are 
perpendicular to the optic axis but they are not 
orthogonal to each other. The y-slope grating is rotated by an angle 0. This rotation reduces the 
amount of y-tilt in the y-slope map by cos(0) and increases the amount of x-tilt in the y-slope map 
by sin(0). Therefore a measurement error is introduced into the Cl and C2 terms in Table 2. 


Cl - ^ sin(0) C2 
* 


006 ( 0 ) 


To determine the effect of such an error, the relationships in Table 2 and Table 3 were 
simulated with a spreadsheet program. The B1 term was set equal to the expression given in 
equation 12 (ideal case), the Cl and C2 terms were set equal to the expressions given in equation 
13 (error case), and all the other B and C terms were set to zero. The misalignment causes an 
error in the amount of Seidel defocus and astigmatism measured. The error in the astigmatism is 
twice the defocus error. The results of the theoretical model for the Seidel astigmatism is shown in 
Fig. 6. Figure 6 compares the magnitude of the error for different grating frequencies and system 
f/#’s. 






Generally, a misalignment is not critical until it alters the amount of aberration measured 
by a quarter of a wave. The sensitivity of the Ronchi test to non-orthogonal gratings is directly 
proportional to the grating frequency and inversely proportional to the system f/#. Figure 7 
compares the experimental results with the theoretical model for non-orthogonal gratings. The 
experimental data agrees with the theoretical prediction. A rotation of ±4* is acceptable for the 
f/2.25 system while a 60 Ip/in grating while only a rotation of ±2* is acceptable for the same system 
with a 100 lp/in grating. 


Effect of Grating Frequency and f/# 

Rotation about the Optic Axis 



— Astigmatism [50 lp/in i JH * 2.25] 

— Astigmatism [1 00 lp/in f/# = 2.25] 

— - Astigmatism [50 lp/in f/# = 4] 


Figure 6 Comparison of different grating 
frequencies and system f/#’s 


Experimental & Theoretical Comparison 
Rotation about the Optic Axis 



x experimen ts Defocut ■ Experimental Astigmatism 
- * Theoretical Defocus — Theor eticaJ Astigmatism 


Figure 7 Comparison of experimental and 
theoretical data 


3.2 Rotation parallel to the grating structure 


To test the effect of rotation parallel to 



Figure 8 Geometry for rotation parallel to 
grating structure 


the grating structure the fixed and variable grating 
are at the focal plane of the system and are 
orthogonal to each other but the variable grating 
is not perpendicular to the optic axis. As the 
grating is tilted the number of fringes in the test 
channel Ronchigram increase. This 
misalignment alters only the Cl term not the C2 
term. The theoretical model for this 
misalignment is derived from the geometry 
shown in Fig. 8. Figure 8 shows the wavefront 
leaving the collection lens and coming to a focus. 
A grating is shown in its aligned position and 
then in its misaligned position. 

In order to determine the number of 
fringes in the Ronchigram, it is necessaiy to 
know the amount of the grating which is 
illuminated. This is determined by calculating 
the length of the line segment ab in Fig. 8. The 
equations for lines 1> 2, and 3 are: 






Line 1 


z-2JI*y 


(14) 


Line! z - -2 Jl* y 


(15) 


Line 3 i - t*n(e) y * « r 


(16) 


The coordinates of points a and b can be found by computing the intersection of the appropriate 
lines. The number of fringes produced in the Ronchigram from the line segment ab is: 


C2 “ “ tanCO) + W * tan(0)F * - t«n(0) W ♦ t®* 0 ) 

Theoretical models were generated with the C2 term equal to the expression given in 
equation 17. The only Seidel aberration that was affected by this misalignment was astigmatism. 
Figure 9 compares the amount of aberration introduced for various grating frequencies and system 
f/#’s. Once again the amount of aberration introduced is directly proportional to the frequency of 
the grating and inversely proportional to the system f/#. However, a change in the system f/# will 
cause a slightly greater change in the amount of aberration introduced than a change in the 
grating frequency would. Figure 10 shows the comparison of the experimental results and the 
appropriate theoretical curve for the system shown in Fig. 4. Both curves have the same form 
although the theoretical curve is slightly higher than the experimental curve. This can be 
accounted for in the inaccuracy of the measurement of the system f/#. The theoretical defocus 
curve lies upon the x-axis. 



Effect of Grating Frequency and f/# 
Rotation Parallel to Grating Struture 



- 60-40 -20 0 20 40 60 

Rotation Angle [degrees] 


- * Astigmatism [50 IpAn V# * 2.25] 

— Astigmatism [1 00 Ip/in t/# E 2.25] 

— Astigmatism [50 tp/in V# * 4] 

Figure 9 Comparison of different grating 
frequencies and system f/#’s 


Experimental & Theoretical Comparison 
Axis Parallel to Grating Structure 



x Experimen tal Pefocm ■ Experimental Astigmatism 
-- Theoretical Defocui — Theoretical Astigmatism 


Figure 10 Comparison of experimental and 
theoretical data 


3.3 Rotation perpendicular to the grating structure 

A theoretical model for the misalignment of the test grating perpendicular to the grating 
structure has not yet been developed. The area of the grating that is illuminated is identical to 
the area of the grating that is illuminated in the case of rotation parallel to the grating structure. 
The difficulty is that rotating the grating in this manner causes the fringes to "fan out", expanding 





Grating Flotation 
Perpendicular to Grating Structure 



Figure 11 Experimental results for rotation 
perpendicular to the grating structure 


on one side of the Ronchigram and compressing 
on the other side. Figure 11 shows the 
experimental results. This misalignment does 
not have a significant effect on the performance 
of the system shown in Fig. 4. 

3.4 Translation along the optic axis 

The fixed and variable gratings are 
perpendicular to the optic axis and are 
orthogonal to each other, but the test grating is 
displaced along the optic axis. Displacing the 
grating towards focus reduces the number of 
fringes that are present in the Ronchigram 
while displacing the grating further away from 
focus increases the number of fringes. This 
misalig nm ent introduces an error into the C2 
term. 


e.v oe.v 

C2 - — ♦ — — (18) 

fl* ft* 

where 6e t ■ the amount of translation along the optic axis. 

This error should affect both the Seidel defocus and astigmatism coefficients. Translating 
the grating further from focus will cause an increase in the amount of astigmatism measured while 
not changing the amount of defocus. This is because the variable channel Ronchigram will have 
more fringes than the fixed channel Ronchigram. The software characterizes the addlt ^ f ^ g6S 
as astigmatism. Translating the grating towards focus will increase the amount of a8 ^ a ^ ls “ 
and decrease the amount of defocus measured. This is because the test channel Ronchigram will 
have fewer fringes than the reference Ronchigram. 

Figure 12 compares the effect of thiB misalignment on systems of various f/#’s and grating 
frequencies. The magnitude of the error is directly proportional to the grating frequencyand 
inversely proportional to the f/# of the optical system. Figure 13 shows the agreement between the 
experimental results for the system shown in Fig. 4 and the theoretical model for that ^ 

translation misalignment along the optic axis becomes important for the f/2.25 system at ±0.5 mm 
for a 50 Ip/in gTating and ±0.25 mm for a 100 lp/in grating. 


3J5 Tr ans lation perpendicular to the grating structure 

Displacing the test grating perpendicular to its grating structure relative to the reference 
grating has no effect on the spacing of the fringes in the Ronchigram^ However * th“ 
displacement can be performed in a controlled periodic motion, for instance # with a piezoelectric 
transducer, then it can be used to implement phase measuring techniques. 


3.6 Tr an slation parallel to the grating structure 

Translation parallel to the grating structure does not have any effect on the spacing or 
orientation of the fringes in the Ronchigrams. 



Effect of Grating Frequency and f /# 
Translation along Optic Axis 



- - Astigmatism [50 Ip/in V# * 4] 

— Astigmatism [50 ip/in f/# * 2.25] 

— A stigmatism [100 Ip/in f/fl * 2.25] 

Figure 12 Comparison of different grating 
frequencies and system f/#’s 


Experimental & Theoretical Comparison 
Translation along Optic Axis 



x I^rtnantaF Drfocu* ■ Eaparrortal Astigmatism 
- * Theoretic al PHoajf — Theoretical Astigmatism 

Figure 13 Comparison of experimental and 
theoretical data 


4. CONCLUSIONS 

One method to fulfill the requirement of two orthogonal slope maps to reconstruct a surface 
„ to J . t w^Lrnel Ronchi inatmment. Tie effect of different greting 

investigated The results presented are based on the use of Zernike polynomia P® 

data and then using linear combinations of those Zernikes to reconstruct the surface shape. Of 
£ £££* SSmento only four led an effect on the 

instrument, the three rotational misalignments and the translation along p 

The severity of these misalignments is dependant upon the frequency of the grating used 
■ fca.t nnd on the f/# of the system. Fast optical systems with higher frequency gratings 
reoufre more c«e in the lament of the gratings thl slower optical systems with low frequency 
gratings The misalignments become critical when they alter the amount of aberration ^asure 
W a afarter of a wave W4. The point at which each misalignment becomes critical depends 
S.„rT^LtoTS-f th. ayrtem and grating being need. The two mort critical 
mtaalignments, independent of system specifications, are translation along the optic axis an 
rotation about the optic axis. 
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Vector formulation for interferogram surface fitting 

David J. Fischer, John T. O’Bryan, Robert Lopez, and H. Philip Stahl 


Interferometry is an optical terting technique that qimntifiea the optical path difference (OK)) between a 
reference wave front and a test wave front baaed on the interference of light. FnngeB are fome ^*“ 
the OPD is an integral multiple of the illuminating wavelength. The resultant tw<>^<^on^pett^ 
is called an interferogram. The function of any interferogram analysis program is to abaci this OPD 
and to produce a representation of the test wave front (or surface). This is accomplished through a 
three- step process of sampling, ordering, and fitting We develop a generalised linear-algebra vector- 
notation model of the interferogram sampling and fitting process 
J&v word * Interferometry, polynomial fitting, optical testing. 


1. Introduction 

Whenever two coherent wave fronts exist at the same 
point in space at the same time, they superimpose. 
The irradiance value at that point is determined by 
the optical path difference (OPD) between the two 
wave fronts. Interferometry uses this phenomenon 
to quantify the shape of a given surface or wave front 
relative to some reference. The two-dimensional 
irradiance pattern created by the OPD between the 
test and reference wave fronts at all pointB is com- 
monly called an interferogram. The challenge is to 
. extract a representation of the test surface (or wave 
front) from such an interferogram (assuming that we 
have perfect knowledge of the reference surface or 
wavefront). This process is accomplished by a three- 
step process: data sampling, ordering, and fitting. 

While the sampling and ordering steps can be 
accomplished by a number of different methods, 
fringe digitization, phase-measuring interferometry, 
or Fourier decomposition, to produce a convenient 
representation of the surface under test, this sampled 
and ordered data (regardless of how they were ac- 
quired) are generally fit to some polynomial set (such 
as the ubiquitous Zernikes). Fringe digitization soft- 
ware has been available to accomplish this process 
since the 1960’s 1 and has been discussed at length for 
Zernike polynomials by using conventional summa- 
tion notation with some matrix notation.^® 

Zer nike polynomials are used because they form a 
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basis set for circular apertures. But they are not a 
basis set for noncircular apertures. In these cases 
other polynomial sets, Buch as Zernike— Tatian 6 or 
Zemike-Mahtgan 78 polynomials (for annulus aper- 
tures) and Legendre® polynomials (for rectangular 
apertures), may provide better results. Ideally an 
interferogram analysis software package should allow 
for easy substitution between these and other polyno- 
mial sets; however, most programs are developed 
exclusively for Zernike polynomials. The first step 
in developing a more general software package capa- 
ble of using any basis set is to develop a general 
Tnn t.hp Tnnt .ical model of the fitting process, completely 
independent of any specific polynomial set. 

In this paper we review the interferogram analysis 
process (sampling, ordering, and fitting) and develop 
a linear-algebra vector notation that clearly and 
concisely describes the interferogram sampling and 
fitting process. The development is completely gen- 
eral (i.e., not limited to Zernike polynomials) and 
assumes only that the polynomials form a basis set. 
The approach performs a least-squares fit of the data 
to the polynomial set, orthogonalizes this result with 
a Gram-Schmidt algorithm, and solves for the final 
weighting coefficients of the initial polynomial set 
with abacksubstitution technique. 

2. Interferogram Analysis 

Because an interferogram has nearly infinite resolu- 
tion, the actual number of points that can be mea- 
sured is limi ted only by the sampling technique. 
There are two main methods for acquiring the data: 
digitization and phase-measuring interferometry. 
With digitization, data points are acquired by sam- 
pling along the center of the light or dark fringes. 
Usually the interferogram is photographed, and the 
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imag e is placed on a digitizing tablet where the data 
points are digitized by a human operator. Alterna- 
tively, the image can be acquired with a video camera 
and digitized by image-processing techniques. Phase- 
measuring interferometry by contrast is an electronic 
process where a video camera samples the wave-front 
phase at every point within its detector array — thus 
creating a high-density, uniform grid of data. An 
illustration of the difference in the two techniques is 
shown in Figs. 1 and 2. In either case the function of 
the sampling process is to acquire data in the form m , 
x,y, where x, y gives the position of the point and m is 
the order of the OPD at that point. 

Once the interferogram has been sampled, the 
acquired data must be ordered. This is because 
interferometry measures the wave-front phase <J> (or 
OPD) in units of modulo 2tr (or modulo X). Thus, 
even though a point in the interferogram represents a 
quantity greater than 2ir, no measured quantity is 
larger than 2tr. For example, an x, y point may have 
a phase of 9ir, but the measured m will equal 0.5, 
which is found by taking the magnitude modulo 2ir. 
Therefore the phase must be artificially recon- 
structed by ordering the data. This is done by 
pncigning integer order number values to each series 
of measured data points. 

For manual or video digitization, data points along 
a given fringe all have the same order, and each 
adjacent fringe has an incrementally larger order 
number. Typically such data consist of integer m 
values, regardless of whether one digitizes bright or 
dark fringes. Assignment of these values is accom- 
plished by the human operator, in the case of manual 
digitization, or the video digitization software itself. 
Alternatively, a phase-measuring interferometer, 
which is not confined to sampling only along a given 
fringe, measures fractional values of m. However, 
its data are still limited by the modulo 2ir condition 
and must be ordered by a process called phase unwrap- 
ping. 

To illustrate the sampling and ordering process, 
consider a one-dimensional parabolic wave front inter- 
fered with a planar wave front to produce an irradi- 
ance pattern (Fig. 3). From the irradiance profile it 
is easy to visualize how an interferogram measures 
the total phase modulo 2 it. Bright fringes occur for 
even multiples of it (i.e., <t> — 0, 2ir, 4i it, . . .), and dark 
fringes occur for odd multiples of -ir (i.e., 4> “ tt, 3ir, 



Fig. I. When an interferogram if sampled by digitization, data 
pointa are usually pla ced along ea ch fringe (either bright or dark). 



Fig. 2. When an interferogram if aampled by phaae-mea*uring 
interferometry, data pointa are placed at each detector location. 
This usually result# in a uniform high -density a a tri pling of the 
pupil. 

6m, . . .). If the pattern is sampled on each dark 
fringe, all the m values will be 0.5. In the typical 
ordering process, these fringes will be assigned values 
of m *= 1, 2, 3, 4, 5. Clearly the result represents a 
parabolic wave front (Fig. 4), but there is a constant 
phase offset between the original wave front and the 
reconstructed data. To extract the original wave 
front, it is necessary to fit the ordered data to a 
parabolic function. Furthermore, if this were a real 
wave front with aberrations, the only way to deter- 
mine their contribution is to perform a general 
polynomial fit. 

Fitting is typically accomplished by a least-squares 
method with a Gram— Schmidt orthogonalization. 
Traditionally interferograms are fit to Zemike polyno- 
mials for several reasons. Ze mikes look like classi- 
cal Seidel aberrations. They provide good results for 
the circular aperture interferogram of most optical 
components. And, since they are orthogonal, they 
are easy to manipulate. Thus it is possible to add or 
subtract coefficients without seriously affecting the 
quality of the fit. In general these coefficients pro- 
vide information about how fiat the surface is, whether 



Horizontal dtotanoa 


Fig. 3. Top curve, One-dim enaional alice of a parabolic wave 
front. Bottom curve, Reaultant irradiance pattern formed when 
♦hi. slice interfere* with a plane wave front. The pluaaes along 
the bottom of the graph represent data pointa aampled at irradj- 
m eg minima. 
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polynomials multiplied by a matrix of Af coefficients: 

f G, l 

[Z r ] ■ [F,(x r . y r ) Fu(x r ,y r )} : • (3) 

[Gy 

Now, if the sampled and ordered data are thought of 
as a sequence of N points, 

rz.i 


HortzortfeJdtfftnoi 


nonwwww* 

Fig. 4, When plotted, the ordered sample point* clearly represent 
m ft a s i hnlir akane. 


the complete set of all N data points can be repre- 
sented as a matrix of polynomials multiplied by the 
coefficient matrix: 


a feature is a bump or hole, whether the part is 
concave or convex, and what aberrations are presen . ^ 

Pi 

3. Vector Formulation w 

o 1 

A. Ideal Surface tl 

The surface height function at any point in the exit 
pupil Z(x, y) can be represented by a linear combina- 
tion of M polynomials F(x, y) and their weighting 
coefficients G: 

Z(x,y)='ZF J (x,y)G J , (D ? 

>-l 8 

where the polynomials are typically a standard set 
' such as the Zemike polynomials. If the coefficients 
are known, the surface function at any pen point is 
found by summing the product of each polynopal ^ 
multiplied by its coefficient. However, m opticd , 

testing the .coefficients are not known. Thus to < 
determine a representation of the si^face function 
one samples aid orders the wave front at many , 

discrete points, fits these data to a polynomial set, and « 

solves numerically for the coefficients. i 

B. Discrete Data 

As discussed in Section 2, the first step in interfero- 
gram analysis is to sample (and order) thecontmuouB 
surface function producing a setof N discrete da 
ooints. With Eq. (1) each sampled point in the pupil 
(x y r ) has a measured value Z t (where r is the sample 
index) that can be represented by a linear combina- 
tion of M polynomials and coefficients: 

Zlx n y T ) - 2 F^y^j, (2) 

where N > M. Alternatively, the value at each 
sample point can be expressed as a matrix of M 
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Fjfo.yi), ... .F^Xi.yillfGi 


i - : : : * (5) 

[Z N \ !/,(**, y N ), .... F*(x n , >n)J L G *J 

To simplify this notation, we define three vectors. 
First, the measured values are defined as a vector 
with N rows. Second, each polynomial Fj eva ^®*^ 
over all N data points is a vector with N rows. Third, 
the M coefficients are a vector with M rows: 

ZA fF;(xi,yi)l G > 

z= : , Fj = : • G- : . • (6) 

Z N \ [f/x n , yw)J l G *. 

With this notation, the sampled data can be repre- 
sented by 

Z - [Fj F*]G. (7) 


C. Fitting to the Polynomials 

The second step in solving for the polynomial coeffi- 
cients is to fit the N data points to the polynomial set 
This is typically accomplished with a least-squares 
Shod, where S is defined as the sum of the square 
of the difference between the data points and the 
fitted polynomials. With the vector notation, S is 

given by 

S*|[Z-[F 1 ,...,F J ,]G]| 2 , ( 8 ) 

or, as more commonly expressed by summation nota- 
tion, 

< 9) 

r - 1 [ jm 1 J 

If a perfect fit were possible, S would be equal tozer^ 
because there would be no difference 
measured values and the representing 
evaluated at the corresponding 
there are always differences between the real surface 
and its representation. Therefore S is nonzero. 


The function of a least-squares fit is to find the 
coefficients for a given set of polynomials, which 

minimizes S. , 

The coefficients are found by taking the derivative 
of S with respect to each coefficient and setting the 
result equal to zero. For example, minimizing S for 
G k gives 

fz.F.fe J,) -it ( 10 ) 

The result of minimizing S for all coefficients can be 
expressed by using summations in matrix form as 


matrix, 


Z-Fjl [F, *~Fi . . . 0 


0 . . . ■ F* 


G, (15) 


and the coefficients can be found by multiplying both 
aideB of the equation by the inverse of the diagonal 
matrix: 

\Lh 1 

II Fi II 2 


XZ r F,(x r ,y r ) 

r-1 


^Z T F^x n y r ) 

r-1 


X *l( X r, yrW\{Xn ?r) ■ ■ • X F*( X r> VrW l(*r, >r) 
r-1 rml 

N N 

X Fl(*n yrWuiXn *) . - • 2 *)**(*» >r) 

_r-l r - 1 J 

G i 

x ; . (ii) 

G*. 

This expression can be greatly simplified with vector 
notation: 

ZF,] rF,-Fi...FrFM* 

: = : . G, (12) 

Z-F*] [Fm-F^.-Fm Yu. 

where the dot product of the vectors F, and F, gives 
the sum of the products of both polynomials over 
every data point: 

F, • F, = X Fix,, y,]Fj(xn *)• ( 13 ) 

r-1 

And the dot product of the vectors Z and F, gives the 
sum of the products evaluated at all data points: 

Z*F ( 14 ) 


However, most orthogonal polynomials are only 
orthogonal over a continuous range, not over discrete 
points. Thus the solution process is not as easy as it 
seems at first. To overcome this limitation, there 
are three possible approaches. The first is to assume 
that the polynomials are orthogonal This may work 
if there is enough sampled data to approximate a 
continuous range, but it can introduce error. The 
second is to take the inverse of the n^dia^omd 
matrix, but this can be numerically difficult and 
inaccurate. Third, the polynomials can be orthogo- 

nalized. T^picaUy, this is accomplished with a uram- 

Schmidt technique. 

4. Grarrv-Schmktt Orthoqonaltzation 

A. Gram-Schmidt Orthogonalization Process 
The Gram-Schmidt orthogonalization process takes 
M arbitrary vectors in an M-dimensional space (M 
space) and creates M orthogonal vectors that form a 
basis for that space. A basis for an M space is 
essentially a coordinate system for that space, bince 
the basis vectors are linearly independent, all other 
vectors in that space may be formed from linear 
combinations of the basis vectors . 

To illustrate the Gram-Schmidt process, we con- 
sider two arbitrary vectors F, and F 2 (Fig. 5), which 
are to be orthogonalized to produce new orthogonal 
vectors 4>. and 4> 2 (Fig. 6). To begin the orthogonal- 
ization process, we define the first orthogonal vector 
to be the first arbitrary vector Fi. Now, with the 
first vector orthogonalized, all other vectors can be 
orthogonalized by projecting the current vector onto 
preceding orthogonalized vector (Fig. 7). Each 


D. Solving for the Coefficients Assuming Orthogonal ^ 
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Fig. 6. Two orthogonal vectors in a two-dimensional space. The 
function of the Gram-Schmidt orthogonalization process is to 
aoate these two vectors. 


original measured surface function values are now 
fitted in terms of the new vectors: 

z - [ 4 *i, .... 4 »*]T. ( 20 ) 

The system is then solved, as shown above, by 
multiplying both sides of the minimized least-squares 
fit equation by the inverse of the resultant diagonal 
matrix. The T coefficients are now known and can 
be used to represent the original surface: 


projection is subtracted from the current vector, 
leaving the orthogonal component. This process is 
continued for the rest of the vectors, each tune 
subtracting off the projections of the current vector 
onto all previous vectors. This creates a set of M, 
orthogonal vectors, which replaces the original set of 

In general, the Gram— Schmidt orthogonalization 
process can be expressed for any vector F, as 


+i m *J 


Cl F 

- vA 


4*. 


.^ 4 >, 4 ». 


4*.. 


(17) 


By letting D Jt , represent the coefficient of 4»„ 


A. 


F, 4>, 
4*. 4», ’ 


(18) 


the orthogonalization of any vector F ; can be ex 
pressed as 


4>. = F j -2d,.4v (19) 

• -1 

These orthogonalized vectors are now substituted 
into the fitting process, replacing the original arbi- 
' trary polynomial vectors. 

Applying the Grarn—Schrnidt Process to the Vector 
Formulation 

From the fitting process, the polynomial vectors, Fj 
through F* [Eq. (7)], are orthogonalized and replaced 
by the new orthogonal vectors, 4>! through 4>*, result- 
ing in the calculation of a new and different set ol 
weighting coefficients. These new coefficients are 
named T. for j in 1-M. Therefore the original coeffi- 
cient vector G is replaced by the vector I\ The 



F5g. 7. Gram-Schmidt orthogonalization process, which sub- 
tracts the projection of the second vector onto the proceeding 
vector to create an orthogonal vector. 


r 


r z • 4 >j 

II 4>! II * 


Z 4>* 

114 ** 11 *. 


( 21 ) 


5. Finding G In Terms of T by Backsubstltutlon 
The polynomials fitted to the measured values are 
typically a standard set, such as the Zemike polynomi- 
als. Therefore, Bince the polynomials do not change, 
the coefficients found in the fitting process of the 
original polynomials can be used to describe the test 
surface and serve as a basis of comparison between 
other test surfaces. However, when the F i through 
F* vectors were orthogonalized, the polynomials r > 
through Fy were modified, creating a new set of 
polynomials, 4>, through 4>*, that are dependent on 
the measured values. Thus the T coefficients that 
compose T cannot serve as a basis of comparison, 
since eveiy Bet of T coefficients is related to a different 
set of polynomials. Because of this, it is necessary to 
find a meanB to convert T to G. 

The first step in finding G in terms of r is to express 
the orthogonalization process by using matrix nota- 
tion. From Eq. (19) the orthogonalization of all M 
vectors can be shown as follows: 


r*i r ' 


F 

T "1 
X 





F 

r 

M 





’ c 

... o' 




D 

21 0 ... 0 

4*i T 


__ 

D 

31 *>32 0 ... 0 





^ • : 

.4 >* T 



[Phi Dh 2 • • • i 



( 22 ) 


Calling the matrix of the D Jt coefficients 2 for simplic- 
ity, we can solve Eq. (22) for the orthogonal vectors, 
■ikm F ift an identity matrix of size M X M . 


4*i r " 


Fi r 

• 

- + 2)' 1 

. 

i 

1 


Jy T . 


( 23 ) 
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Recalling Eqs. (7) and (20), we can dearly see that two 
equivalent representations of the surface exist, one 
using the original polynomial vectors and the other 
using the orthogonalized polynomial vectors. There- 
fore the following statement can be made: 

[Fj Fjf]G * [4>j 4V]T. (24) 

Now, if we rewrite Eq. (23) using transposes 

f*i.."A»F-t^ + *)' , lF| l F^f, (25) 

and the transpose of both sides of the resultant 
equation is taken, an expression for the matrix of 4>j 
through in terms of F, through F* is produced: 

{♦i, . . . ,«*] « [Fj, .... T M l(S + 2)-') t . (26) 

Substituting Eq. (26) into Eq. (24), 

[F if ...,F*]G ■ [Fj F 'mI(S + 9)-'YT, (27) 

and dividing both sides by the matrix of the polyno- 
mial vectors provide the desired solution of G in 
terms of the D ^ coefficients and T : 

G - [(-f + (28) 

6. Conclusion 

Interferometry is an optical testing technique in 
which information about a test surface is extracted 
from an interferogram. This is accomplished with a 
three-step process. First, the data are sampled, usu- 
ally by one of two main methods: digitization or 
phase-measuring interferometry. Second, the sam- 
pled data points are ordered. Third, the ordered 
data are fit to a polynomial set. Typically, Zemike 
polynomials are used because they form a basis set for 
circular apertures. However, they are not a basis set 
for noncircular apertures. In these cases, other 
polynomial sets such as Zemike-Mahqjan (for annu- 
lus apertures) and Legendre (for rectangular aper- 
tures) may provide better results. In this paper we 
developed a completely general linear-algebra vector- 
notation model to describe the polynomial fitting 
process independent of the polynomial set. The 
basic process is to describe the sampled data as a 
linear vector: Fit that data to an arbitrary polyno- 
mial set by using a least-squares method. Orthogo- 


nalize the polynomials to the data vector by using a 
Gram-Schmidt technique. Solve for the original 
polynomials weighting coefficients by using a backsub- 
etitution technique. 

From an implementation point of view, by treating 
the polynomial set as a vector, it is possible to write a 
algorithm that hf»* no knowledge of the actual polyno- 
mial set being fitted. Thus it is simple to change 
between polynomial sets for various aperture shapes. 
Also, it is not necessary to actually solve for the 
orthogonal polynomial coefficients to arrive at the 
original polynomial weighting coefficients by using 
the backsubstitute technique. 
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ABSTRACT 

One approach for obtaining a surface representation is to fit Zeraike polynomials (in a least squares sense) to 
discrete data points in the full aperture. The mathematics for this have been described using both matrix 1 and vector 
notation 2 . Additionally, vector notation has been used to describe bow to obtain a surface representation from orthogonal (x, 
y) slope data 3 . The result of that paper was a matrix operator for linearly combining the first eight Zernike polynomial 
coefficients fit to x- and y-slope data to produce a Zernike polynomial surface representation. This paper extends that 
process by presenting a systematic approach for obtaining the linear relationship between slope and surface using the first 49 
Zernike polynomials. 


1. INTRODUCTION 

Full-aperture interferometry techniques are used to measure surface characteristics of a test object via non- 
destructive optical methods. 1 Traditionally, Zernike polynomials are fit (in a least squares sense) to the measured data to 
produce a standardized representation of the surface. The parameters which describe the surface’s characteristics are, the 
weighting coefficient of the Zernike polynomials. This fitting process has been explained in detail using both matrix and 
vector analysis 1,2 . 

While surface interferometry is be simple to analyze, it is not always used; in some cases it is better to measure 
slope, such as the Ronchi test does. Because, these techniques yield slope data in a particular direction, to create a surface 
representation, it is necessary to measure slope in two orthogonal directions: in the x- and y-directions. The data is not fit to 
the derivatives of the Zernike polynomials because the derivatives are not orthogonal. Rather, to reconstruct the surface, 
these orthogonal slope data sets are fit to Zernike polynomials to create a representation of the test object's slope, and 
combined according to a linear operator. 3 This paper describes the process of creating a polynomial surface representation 
using the first 49 Zernike polynomials based on a least squares fit of Zernike polynomials to orthogonal slope data. 

2. FITTING ZERNIKE POLYNOMIALS TO THEIR DERIVATIVES 

The first step in creating a Zernike polynomial, surface representation based on fits of orthogonal slope data is to 
determine which linear combinations of Zernike polynomials will recreate the x- and y-derivatives of the Zernike 
polynomials. To do this, it is first necessary to examine the Zernike polynomials. 

2.1 Description of Zernike polynomials in R terms 

The Zernike polynomials are a basis set, orthogonal over a unit circle. There is an infinite number of these 
polynomials, since they are constructed from a generation formula, but in many applications only die first 36 or 49 
polynomials are used They are commonly used for fitting data from interferometric applications because they represent the 
Seidel aberrations well. 4 In this paper, two different notations will be used to represent the Zernike polynomials. They are 
the Z. notation and the Rf” notation. In the case of the Z, notation, die subscript, n, indicates which polynomial is 
lying represented The subscript is a positive, integer value starting at zero. Thus, Z, 0 refers to the eleventh Zernike 
polynomial. For the Rf” notation, r is the radial magnitude of the polynomial and m is the radial frequency. The sign of m 
also designates whether the trigonometric term involved in the Zernike polynomial is a sine or a cosine: m £ 0 designates a 
cosine and m < 0 designates a sine. Finally, r £ 0 always, andr^lml always. These cases are summarized below: 


£?££ . 


Vet. <2.°o3 




/ J%Z. 




In the remainder of this paper, the Zemikes will be considered using rectangular coordinates (x and y), not polar (r and 6). 
Examples of this notation follow: 


R 4 2 = Z, 2 = (4^ -3) r 2 sin (20) = 8x 3 y + 8xy 3 - 6xy 
Rj = Z, = r 3 cos (30) = x 3 - 3xy 2 
Figure 2. Illustration of notation of Ze mike polynomials 


23 Selection of terms to fit 

The next step is to create the x- and y-derivatives from the Zemike polynomials. The data is not directly fit to the 
derivatives of the Zemike polynomials because the derivatives are not orthogonal. For each of the first 49 Zemike 
polynomials, use linear combinations of Zemike polynomials. Instead of trying a combination of all 49 polynomials to 
determine which polynomials construct a given derivative, the Rf" notation dearly indicates which polynomials to include 
in the linear combination. 


The first step is to determine the highest order term in the summation. This is found according to the following 
rules given in Figure 3, below: 


3R r B 

Tx 

use R“ V 

unless m+1 > r-1, then use 

nlD-l 

K r-1 

3R r m 

5y 

MCA + 

use K r _ j 

unless m+1 > r-1, then use 

p-(m-l) 

K r-1 


(a) m £ 0 


3R" 

Tx 

use R“ V 


nB4] 

K r-i 

unless hn-ll > r-1, then use 

3R" 

Ty 

use R^*” 

unless hn-ll > r-1, then use 

K r-i 


(b) m < 0 


Figure 3. The four cases for determining the highest order term of a linear summation 
which produces the derivative of the Zemike polynomial specified by the initial R term. 




The rest of the terms used in the linear summation are found from the highest order term. This is shown below, with 
an example following. 


3R“ 

To find 

for m £ 0,even, use 

nin-2 pn-4 pO pm pm-2 p0 pO 

» ■ • •? •'►j f _ 2 * — 2 * *•'» — 2* "** 

9R“ 

To find ^ 

for m > 0, odd, use 

pm-2 pm -4 pi pm pm-2 pi pi 

* R r * * v r t *S-2* *S-2 * 1v t-2* **** **1 

9R“ 

To find ^ 

for m < 0,even, use 

pm + 2 pm +4 p 2 pm pm — 2 p-2 

t i "•» ’ ‘S 1 — 2* *S'-2 * — 2* •“» * v o 

9R“ 

To find ^ 

for m < 0, odd, use 

pm + 2 pm+4 p-1 pm p» + 2 p-1 p-1 

X r » » •••! I'j t •S 1 — 2* *^-2 * •••* — 2* *••• 

Figure 4. Selection of all terms used in a linear summation to produce a given derivative. The initial R 
value is the one chosen according to Figure 3. It should be recalled that m £ r always, so R terms involving 
a superscript value greater than its subscript will not exist, and are not used in the linear summation. 


To find the x-derivative: To find the y-derivative: 


(a) 

Z n = 

Z l3 = Rj 

(b) 

R 1^*1 

r’-*r; 2 

(c) 

p2 pO p2 pO pO 
IV 2 t *^0 

** 
* i 

KJ 

w 

(d) 

Z,p Z 8 , Z 4 , Zy Z 0 

Z 12 , Zj 

(e) 

dRj 

— bjiZjj + bgZ g + b 4 Z 4 + b 3 Z 3 + boZo 

3Rj 

- CjjZjj + CjZj 


Figure 5. An example of selecting the terms of the linear summation, (a) The Zemike 
polynomial being considered, (b) Finding the highest order term of the summation for 
the x- and y-derivatives. (c) All terms involved in the linear summation for the creation 
of the x- and y-derivatives. (d) The terms expressed in Z„ notation, (e) The linear 
summation used to create the derivatives; bn, bg, b 4 , bo and Cj2 and c$ are constants. 

This process can be repeated for the x- and y-derivatives of the first 49 Zemike polynomials. 

23 The fitting process 

With the polynomials necessary to create each derivative now chosen, it is possible to find the amount of each 
polynomial needed. F-arh polynomial in the linear summation is multiplied by an arbitrary constant, as illustrated in Figure 
S. To give the exact representation of the derivative in terms of Zemike polynomials, numerical values for these constants 
were found using a computer algebra application, MAPLE V, on a NeXT computer running NeXTSTEP v3.0. Two 
approaches to the fitting were used: for simpler summations a simple MAPLE procedure worked well, but for certain higher 
order terms more effort is involved. 

The simple MAPLE procedure begins with a linear sum of the Zemike polynomials that will create the derivative. 
Each Zemike in the sum is multiplied by a constant. Then, all coefficients are matched between the linear sum and the 
derivative. (Coefficient matching is the process of finding the coefficients of a given power of x or y on both sides of an 





equation and setting these terms to be equal.) The set of equations generated by coefficient matching art solved for the 
constants of the summing polynomials. The MAPLE V codes for doing this for the x- and y -derivatives of the 13th Zemike 
ait shown below. 


Consider the 13th Zemike 


The linear sum of 
Zemike polynomials 


sum :* bO*zO + bl # z3 + b2*z4 + b3*z8 + b4*zll: 
sum : * expand (sum) : 

d : *13 : ^ 

linequset :■{} : coeffset :* {}: 
for k from 0 to 4 do 

•create equation matching x-coeffs of x-derivative to x-coeffs of sum 

linequx.k :* coeff (zx.d, x, k) « coeff (sum, x, k) ; ^ “ 

•create equation matching y-coeffs of x-derivative to y-coeffs of sum 
linequy.k :* coeff (zx.d, y, k) ■ coeff (sum, y, k) 

•add the equations to the set of equations 
linequset := linequset union {linequx.k, linequy.k}? 

•add the coeff to the set of coeffs 
coeffset coeffset union (b.k); 

od: 

# solve for the coeffs 
solns : =solve ( linequset, coeffset) : 


Solve the equations 
for the weighting 
coefficients. 


solns : = {b4 = 5., b3 = 5., b2 


bl = 3., bO = 1.) 


Match the 
-coefficients to 
powers of x 

Match the 
^coefficients to 
powers of x 


The values of the 

weighting 

coefficients 


Figure 6. The MAPLE procedure used to find the weighting coefficients of the linear sum to create the x- 
derivative of the 13th Zemike polynomial. The coefficients b4 .. bO correlate to b^, bg, b 4 and b^ respectively 


sum := c0*z5 + cl*zl2: 
sum : = expand (sum): 

d:=13: 

linequset :*{} : coeffset 
for k from 0 to 3 do 


Consider the 13th Zemike 


The linear sum of 
Zemike polynomials 


:■ U 


•create equation matching x-coeffs of x 
linequx.k coeff (zy.d, x, k) * coeff (sum, x. 


derivative to x-coeffs of sum 
k); 


Match the 
- coefficients to 


•create equation matching y-coeffs of x-derivative to y-coeffs of sum powers of x 
linequy.k :* coeff (zy.d, y, k) * coeff (sum, y, k) 

•add the equations to the set of equations ^ Match the 

linequset : = linequset union (linequx.k, linequy.k); ‘^^-^.coefficients to 

•add the coeff to the set of coeffs powers of x 

coeffset := coeffset union (c.k); 
od: 

• solve for the coeffs 
solns : *solve (linequset , coeffset ) 


Solve the equations 
for the weighting 
coefficients. 


solns := {cl * 5., cO = 3.} 


The values of the 

weighting 

coefficients 


Figure 7. The MAPLE procedure used to find the weighting coefficients of the linear sum to create the y- 
derivative of the 13th Zemike polynomial. The coefficients cl and cO correlate to Cj2 and c 5 , respectively. 


Unfortunately, for some cases, this process does not work. The above procedure does not produce integer values for 
the coefficients. When this happens, it is necessary to use more finesse. 


The procedure starts by first identifying the highest order terms in the derivative. The coefficients of these terms are 
matched to the equivalent terms of the highest order summing polynomials, and the weighting coefficients are solved for. 
The highest order polynomials are multiplied by the weighting coefficients, and subtracted from the derivative. The process 
is continued by matching the highest order terms of the remainder to the highest order terms of the remaining summing 
polynomials and solving for the corresponding coefficients. By iterating this process, all weighting coefficients will be 



found. This is illustrated for the x-derivative of the 41st Zeraike: 

> zx41; 

6 43 5 2 7 4 5 

784.0 y x + 560.0 x y - 336.0 y x - 112.0 y - 840.0 y x ♦ 168.0 y 

2 3 

+ 180 . 0 y x - 60 . 0 y 

> expand (z3 9 ) ; 

6 43 52 7 4 32 5 

35.0 y x - 35.0 x y - 63.0 y x + 7.0 y - 30.0 y x ♦ 60.0 y x - 6.0 y 

> expand ( z30) ; 

6 43 52 7 4 32 

63.0 y x + 105.0 x y + 21.0 y x - 21.0 y - 90.0 y x - 60.0 y x 

5 2 3 

+ 30.0 y + 30.0 y x - 10.0 y 

> expand (z23 ) ; 

6 43 52 7 4 32 

35.0 y x + 105.0 x y + 105.0 y x + 35.0 y - 60.0 y x - 120.0 y x 

5 2 3 

- 60.0 y + 30.0 y x + 30.0 y - 4.0 y 

> expand { z26) ; 

4 3 2 5 

5yx - 10 y x + y 

> expand ( zl 9 ) ; 

4 3 2 5 2 3 

15.0 y x + 10.0 y x - 5.0 y - 12.0 y x ♦ 4.0 y 

> expand ( zl4) ; 

4 3 2 5 2 3 

10.0 y x + 20.0 y x ♦ 10.0 y - 12.0 y x - 12.0 y + 3.0 y 

> expand ( zlO) ; 

2 3 

3 y x - y 

> expand ( zlO); 

2 3 

3 y x - y 


> expand ( z7 ) ; 



> expand (z2 ) ; 


2 3 

3 . 0 y x + 3,0 y - 2,0 y 


y 

expand (z2 ) ; 

It is seen that the highest order terms in the derivative of Z41 are of magnitude seven: yx 6 , x 4 y * , yV* y 7 and there 
are also three summing polynomials with terms of order seven: Z39* Zjq, Z23. Thus* the first three weighting coefficients are 
found by using the coefficients of the yx 6 , x 4 y* t and y*x 2 terms: 

> solve ( (35*b8+63*b7+35+b6 • 784, -35*b8*105*b7*105*b6 • 560, -63*b8+21*b7+105*b6-- 
336), 

(b8, b7, b6) ) ; 

> 

(b7 . 8, b6 - 0, b8 - 8) 

> aaaignC) ; 

> stuff t * «xpand(sx41 - (b8+s39 ♦ b7*x30 ♦ b6*s23))j 

4 5 2 3 

•tuff t- 120.0 y x - 24.0 y - 60.0 y x ♦ 20.0 y 
With the first three weighting coefficients found* the process is continued for the yx 4 , y^x 2 , and y 5 terms. 

> •olv#( (5*b5+15*b4*10*b3 - 120, -10*b5+10*b4+20*b3 ■ 0, l*b5-5*b4*10*b3»-24) , <b5, 
b4, b3 ) ) ; 

> 

(b3 »0, b4»6, b5*6) 

> aaalgn(") ; 

> 0tuff2 1 ■ axpand ( atuf f * (b5+z26 ♦ b4*xl9 ♦ b3*zl4)); 

2 3 

atuff2 im 12.0 y x - 4.0 y 

The last of the weighting coefficients are found by considering the yx 2 * and x 3 terms. 

> solve ( (3*b2+3+bl«12, -l*b2+3*bl— 4) , (bl, b2>); 

(bl - 0, b2 - 4) 

> azsignC) ; 

> »tuff3 t • expand ( a tuf f2 - (b2*zl0+bl*x7) ) ; 

•tuff 3 t« 0 

The final coefficient, aO, must be equal to zero. 

Using the described methods, the constants are found for all polynomials in all linear combinations for the x- and 
y-derivatives of the first 49 Zemike polynomials. 

3 . INITIAL APPROXIMATION OF THE SURFACE FROM THE DERIVATIVE FIT 

When the orthogonal* slope data is (te&st squares) fit to Zemike polynomials, Zq ... Z^, a set of weighting 
coefficients for the x-slope and the y-slope are produced. These weighting coefficients will be termed: B 0 ... B43 for the x- 
slope, and Cq ... C4& for the y-slope. The weighting coefficients are now used to generate a surface representation composed 
of Zemike polynomials. 



To find the weighting coefficient for a given Zeraike polynomial, first find the largest constants) used to recreate its 
x- and y -derivatives. The coefficients from the least squares fit corresponding to the largest summing constants should then 
be added together, and divided by the sum of the largest constants. This gives an approximation to the weighting coefficient 
of the Zemike polynomial to represent the surface. 

The 13th Zemike polynomial, Z| 3 , is considered for an example. The largest constants of the x-summatioo are bg 
and b|], with a value of 5. (See Figure 6) The polynomials corresponding to these constants are Z \ ] and Zg. Likewise, for the 
y-derivative, that the largest constant is Cj 2 . also with a value of 5, corresponding to Z 12 . To reconstruct, add the x-slope 
weighting coefficients (Bu and Bg) and the y- si ope weighting coefficient (C| 2 >. This sum should then be divided by 
a3+a4+a2 = 15 to normalize. Thus, the initial approximation of the weighting coefficient of Zj 3 based on the slope data is 
Aj 3 = [Bu + Bg + C 12 I/I 5 . 


4. IMPROVED APPROXIMATION OF THE SURFACE FROM THE DERIVATIVE FIT 

The approximation of the Zemike polynomial representation of the surface found above is improved by subtracting 
off higher order terms because the higher order terms include lower order terms. The higher order terms which should be 
removed are those of the same radial frequencyof, but with greater radial magnitude than, the Zemike polynomial being 
considered. These terms are found from the R, notation, as shown below: 


(a) For a given Zeraike polynomial Rf" , remove the coefficients of: r“ + 4> r“ + 4> . . . 

(b) For the Zemike polynomial Rj , remove the coefficients of: R^R^Rj, 

For the Zemike polynomial Z 13 , remove the coefficients: A a , A J3l A^ 

Figure 8. (a) Determining higher order terms to subtract from surface approximation. 

(b) An illustration of this procedure using Z t3 . 


The resulting weighting coefficients are thus: 



AO 

= 

0 


A1 

St 

[BO] - A6 - A13 - A22 - A46 

A2 

= 

[CO] -hi - A14 - A23 - A3 4 

- A47 

A3 

E 

[B1 + B2]/8 - A8 - A14 - A23 - A34 

A4 

= 

[B1 - Cl]/4 - All - A20 - A31 - A44 

A5 

St 

[B2 + Cl]/4 - A12 - A32 - A45 

A6 

= 

[B4 + B3 + C5J/9 - A14 - A23 - A33 - A46 

A7 

E 

[B5 - C4 + C3]/9 - A14 - A23 - A34 - A47 

A8 

E 

[B6 + C7J/16 - A15 - A24 - 

A3 5 - A48 

A9 

E 

[B4 - C5] / 1 6 - A18 - A29 - 

A42 

AlO 

E 

[B5 + C4]/16 - A19 - A30 - 

A43 

All 

s 

[B9 + B6 - C7 + C10]/16 - A20 - A31 - A44 

A12 

E 

[BIO + B7 + C6 - C9]/16 - A21 - A32 - A45 

A13 

= 

[Bll + B8 + C12]/15 - A22 - 

A33 - A46 

A14 

E 

[B12 - Cll + C8]/15 - A23 - 

A3 4 - A47 

A15 

JC 

[B13 + C14J/24 - A24 - A35 

- A48 

A16 

E 

[B9 - C10J/8 - A27 - A40 


A17 

= 

[BIO - C9]/8 - A28 - A41 


A18 

= 

[B16 + Bll - C12 + C171/20 

- A29 - A42 

A19 

= 

[B12 + B17 - C16 + Cll]/20 

- A30 - A43 

A20 

B 

[B18 +B13 - C14 +C19]/24 - 

A31 - A44 

A21 

= 

[B19 + B14 - C18 + C13]/24 

- A32 - A44 

A22 

S 

[B20 + B15 + C21J/21 - A34 

- A46 

A23 

s 

[B21 + C15 - C20]/21 - A35 

- A47 

A24 

= 

[B22 + C23J/32 - A35 - A48 


A25 

= 

[B16 - C17]/10 - A3 8 


A26 

B 

[B17 - C16J/10 - A3 9 


A27 

E 

[B25 + B18 - C19 + C26]/24 

- A40 

A28 

= 

[B19 + B26 - C25 + C18J/24 

- A41 

A29 

— 

[ B27 + B20 + C28 - C21J/28 

- A42 

A30 

E 

[B28 + B21 - C27 + C20J/28 

- A43 

A31 

E 

[B29 + B22 - C23 + C30]/32 

- A44 

A3 2 

= 

[B30 + B23 - C29 + C22]/32 

- A45 

A3 3 

S 

[B24 + B31 + C32]/27 - A46 


A3 4 

= 

[B32 + C24 - C311/27 - A47 


A3 5 

E 

[B33 + C341/40 


A3 6 

E 

[B25 - C26J/12 


A3 7 

s 

[B26 + C25]/12 


A3 8 

S 

[B36 + B25 + B37 - C28]/28 


A3 9 

E 

[B37 + B27 + C27 - C36]/28 


A40 

s 

[B38 + B29 + C39 - C30]/32 


A41 

E 

[B39 + B30 - C38 + C29]/32 


A42 

E 

[B40 + B31 + C41 - C32]/36 


A43 

S 

[B41 + B32 - C40 + C31]/36 


A44 

S 

[B42 + B33 + C43 - C34]/40 


A45 

E 

[B43 + B34 - C42 + C33]/40 


A46 

E 

[B44 + B35 - C45]/33 


A47 

E 

[B45 - C44 + C35]/33 


A48 

E 

[B46 + C47]/48 


Figure 9. Weighting coefficients to first 48 Zeroike polynomials to reconstruct a 
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5. CONCLUSION 


The means to obtain a surface representation of by fitting Zernike polynomials (in a least squares sense) to discrete 
data points in the full aperture has been described in the past via matrix and vector notation. Additionally, vector notation 
has been used to describe how to obtain a surface representation from orthogonal (x, y) slope data 3 . The result of that paper 
was a matrix operator for linearly combining the first eight Zemike polynomial coefficients fit to x- and y-slope data to 
produce a Zemike polynomial surface representation. This paper extended that process by presenting a systematic approach 
for obtaining the linear relationship between slope and surface using the first 49 Zemike polynomials. 
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ABSTRACT 

Ronchi interferometry is an optical testing technique similar to shearing interferometry. A coherent wavefront is 
interfered with a sheared form of itself by placing a periodic grating at, or near, the focus of an optical system. The resultant 
interference pattern contains information about the wavefront’s slope in a direction perpendicular to the grating structure. The 
wavefront can be reconstructed from two orthogonal slope data sets via the process of sampling, ordering and fitting. This 
paper develops a linear-algebra vector notation model of the interferogram sampling and fitting process. 

1. INTRODUCTION 


Whenever two coherent wavefronts exist at the same point in space, they superimpose, to create an interference 
pattern. The irradiance at any point, E(x,y), as determined by the Optical Path Difference (OPD) between the two wavefronts is 
described by: 

E(x,y) = Eq 1 + cos (^OPD)J ( 1 ) 


\ = measurement wavelength 

( 2 ) 

The two dimensional irradiance pattern created by the optical path difference at all points is called an irradiance pattern. 

For a sheared interferogram, the OPD is given by the difference between the wavefront and its sheared self. This difference is 
called the differential wavefront. 

OPD = W (x, y) - W (x + Ax, y) = AW (x, y) (3) 

By rewriting the differential wavefront as: 

8W (x, y) 

AW(x, y) = ^ Ax = a(x, y)Ax ( 4 ) 


the wavefront slope can then be written as: 

a(x, y) = 

Therefore, the irradiance pattern can be expressed as: 


8W(x,y) AW (x, y) 


E (x, y) = E 0 1 + cos (^< 

L 


X = — = equivalent wavelength 
cq Ax ° 


The Ronchi test is similar to shearing interferometry, only the shear is caused by a periodic grating placed at, or near, 
the focus of the optical system under test. Like shearing interferometry, the resultant irradiance pattern (Ronchigram) contains 
slope information in a direction perpendicular to the shear direction - the direction of the grating structure. From this 
similarity, it can be shown that the Ronchigram can be expressed by: 

E(x,y) = E 0 l+cos(^a)J (8) 


SP2E t'c/, 17 SS , /??<*, 



where: 


a = wavefront slope 

x - d 
*« 2F/# 


(9) 

( 10 ) 


( 11 ) 

( 12 ) 


d = Grating Spacing 
F/# = Test Beam F/# 

Thus, by proper scaling of the interferogram, conventional interferogram analysis techniques will provide a direct 
measurement of the wavefront slope. The challenge is to extract a representation of the lest surface from such an 
interferogram. This will be done using a linear- algebra vector notation model of the interferogram sampling and fitting process 

2. INTERFEROGRAM ANALYSIS 

fiuino ^ onve " ,io " a ' wavefron, . in,erfer °g ram ^alysis is accomplished via a three step process: data sampling, ordering and 
The result of his process is a representation of the wavefront under lest. For a shear interferogram, such as a 
Ronch.gram, the analysis provides a representation of the wavefront’s slope in a particular direction (such as the x- or y- 

slope representations 3 repreSema "° n ° f ,he ac,ual wavefront from a Ronchigram, it is necessary to combine two orthogonal 
2.1 Sampling 

I . ® ecause an interferogram has nearly infinite resolution, the actual number of points that can be measured is limited 

I rfl e a T I 8 '' ChniqUe , T Cre ^ tW ° main meth0dS f ° r acc l uirin g the d ata: digitization, and phase-measuring 
, h l f ? r ^ ^ dlg, ' lzatl0ni data P° ,nls are acquired by sampling along the center of either the light or dark fringes Usually 
the interferogram is photographed and the image is placed onto a digitizing tablet where the data points are dioitized by a 
human operator. Alternatively, the image can be acquired with a video camera and digitized using imaoe-processino 

Dhtenr 5 ' meaS r ng m ; erferometr >’ b >' c °ntrast, is an electronic process where a video camemsamples the wavefront 
phase every point within its detector array - thus creating a high density, uniform grid of data. An illustration to show the 
difference in the two techniques is shown in Figure 1 and Figure 2: 



Fig. I An illustration of the placement of data points 
w hen an interferogram is sampled by digitization 



Fig. 2. An illustration of the placement of data points when an 
interferogram is sampled by Phase-Measuring Interferometry 



For this paper, it is assumed that the data is sampled by manual digitization along the dark fringes. The data is 
recorded in the form (m,x,y) where (x,y) gives the position of the point and (m) is the order of the sampled fringe - which as 
will be discussed in the next paragraph, is related to the height of the surface in units of 2n. 

2.2 Ordering 

Once the data has been sampled, it must be ordered. Because the irradiance pattern is proportional to cos(k*OPD) the 
data is in units of modulo 2 n or modulo Thus, even though a point in the interferogram represents a quantity greater than 2k, 
no measured quantity will be larger than 2n. For instance, an (x,y) point may have an (m) that is of magnitude 9k. The 
measured (m) will equal 7t, which is found by taking the magnitude modulo 2n. Therefore, the magnitude of the data must be 
artificially reconstructed by ordering the data. This is done by assigning integer order number values to each series of 
measured data points. For example, all data points along a given fringe will have the same order and each adjacent fringe has 
an incrementally larger order number. 

For manually digitized data, this ordering is accomplished by the trained operator. For video digitization, this is 

accomplished by image processing techniques. And for phase-measuring, this ordering is accomplished by a process called 
phase unwrapping. 

2.3 Fitting 

Once the data has been sampled and ordered, it must be fitted to a polynomial set to represent the shape of the test 
surface, generally, this is accomplished using a least squares method with a Gram-Schmidt orthogonalization. From this fit 
several questions can be answered: how fiat is the surface, are the features concave or convex, and what aberrations are present. 

Traditionally, Zemike polynomials have been used for this representation for several reasons. They look like classical 
Seidel Aberrations. They provide a good fit to the circular aperture interferogram of most optica! components. And, because 
they are orthogonal, they are easy to manipulate. It is possible to add or subtract coefficients without seriously affecting the 
quality of the fit. c 


3. VECTOR FORMULATION 


3.1 Purpose of the vector formulation 

The purpose of this section is to summarize simply and clearly the analysis process of a Ronchigram using 
general linear-algebra vector notation, and to develop a generalized mathematical model with which to explore the effects of: 
sampling parameters, fitting polynomials, and the surface-fitting algorithm on the accuracy of the reconstruction of the test 
surface. The remainder of this paper presents a general vector based formulation for the surface-fitting process. 

3.2 Representing the measured values 

The wavefront function at any point in the exit pupil, W(x,y), can be represented by a linear combination of M 
polynomials, Fj(x,y), and their weighting coefficients, Gj, for j in 1 to M. 

M 

W(x,y) = ^G^fx.y) (13) 

j- i 

Similarly, the wavefront slope at any point in the exit pupil, a(x,y), can be represented by: 

M 

a(x,y) = ^GjVFjfx.y) (14) 

j« 1 

Unfortunately, the Ronchi test does not make general slope measurements. It makes independent measurements of x-slope and 
y-slope which can be represented as: 



a 


(15) 


M 

,(*.y> = X G *i (F *i (x ’ y)) 


where: 


M 

a y (x, y) = X G yj (F yi (x,y)) 

j- 1 


F xj (x,y) 
Fyj(X. y) 


5Fj(x,y) 

6 x 

SFjtx.y) 

5 y 


It is therefore necessary to perform two separate fitting operations to determine G X j and G y j . 


(16) 


(17) 

(18) 


A sampled and ordered x-slope interferogram data set be thought of as a series of N measured data points, where N > 
M. At any point in the interferogram, (x p y r ), there is a measured value, o^p where r is in 1 to N, the number of sampled points. 
The r Ui measured x-slope value is then written as: 

M 

«„(x f .y r ) = X G *J (F *j (x,y)) where N > M ( 19 ) 

j= i 

Alternatively, the single measured value, a XP can also be expressed as a matrix of polynomials multiplied by a matrix of 
coefficients. 


[ a J = [F„ ( X r-y r ) ■ F »M( X r-y r )] 



(20) 


Thus, the entire series of N, measured x-slope values can be represented as a matrix of polynomials multiplied by the 
coefficients matrix. (See Equation 21) 


a 


Xl 


a 


xNl 


F«i (x,.y,) ... F xM (x, t y,) 

I I 

F xl (X N , y N ) ... F xM (x N .yN) 


J xl 




( 21 ) 


In order to simplify this notation, the measured values can be considered to be a vector with N rows. Likewise, the 
coefficients can be viewed as a vector with M rows, and a given polynomial evaluated at all N data points is a vector with N 
rows. 




o,r 

J 

F »j(X|.yi) 

1 

«XN 

G, = 

i 

G xM 


i 

. f ,j(x n . y N )_ 


Using this notation, the fitted values can be represented by: 



( 22 ) 


( 23 ) 


3.3 Finding the best fit 

The polynomials can be fit using a least squares method, where S is defined as the sum of the square of the difference 
between the data points and the fitted polynomials. If a perfect fit were possible, then S would be equal to zero, because there 
would be no difference between the measured values and the representing polynomials evaluated at the corresponding points. 





b " w “?. lhe •*.-** «- 1“ Therefore, S win be ttotvrere, The 

ynomials, which minimizes S. 
cpressed using summations. (S' 

Hfc~[Kj ... <g«.]| 


. ' w..^ uic icoj bunace ana ns representation I here! 

hrncMon ofa leas, squares fit is lo find the coefficients, fo, a given se, of polynomials, which minimiees S. This process can be 
described nsrng vector no, anon (See Equation 24), as more commonly expressed using summations. (See Equation 25? 


(24) 


M 


S= X a *r-X G »j F *j (X r’yr> 

r- I L j- 1 


(25) 


The coefficients are found by taking the derivative of S with respect to each coefficient and setting the result equal 
example, mini mizine S for C . 6 tltauuc M u ^ 


zero. For example, minimizing S for G xk gives: 

N 


to 


N M 


= XX G *i F *i (x '’yr)F* k (x r .y r ) 

r * 1 r * 1 j * 1 

The result of minimizing S for all coefficients can be expressed using summations in matrix form as: 

N N 

X F ..(Vy,)F„(X,y r ) •• X F *M (X r’ y r )F *l (X r’yr) 
r * 1 r= 1 

I I 

N N 

X F *' (X r’yr) F *M(X r .y,) ... X F -(Vy r )F xM (X,y r ) 
J ■= I r. 1 


(26) 


X a *r F *t( X r> y,) 


i= I 


I 


X“*r F *M( X r. y r ) 


Lr= 1 


J 




L G xmJ 


(27) 


This expression can be greatly simplified with vector notation (See Equation 28), where the dot product of the vectors Fxi and 
*j gives the sum of the products of both polynomials over every data point. (See Equation 29) 


«,*F X1 

| 

II 

F>i • F xl . 
| 

•• F xl . F xM 

i 



1 

jx a 

1 

* FjM # FxM 

a x • f xM _ 


F*m • F xj . 


F «i • F xj = £ F xi (x r , y r ) F xj (x f , y r ) 

r * ] - 

Likewise the dot product of the vectors « x and F xj gives the sum of the products evaluated at all data points. 

N 

“x ,F ‘j = X a *. F *j (X r-yr) 


(28) 


(29) 


(30) 


r = 1 


3.4 Solving for the coefficients assuming orthogonal polynomial vectors 

If the F vectors are orthogonal over the sampled data points, then the minimization process produces a diagonal 
matrix (See Equation 31), and the coefficients can be found by multiplying both sides of % equation by the inverse of the 
diagonal matrix. (See Equation 32) The process is repeated using the y-slope data to find Gy 
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r- 
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a « * F* i 


Fxi • F x i . 
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Ji 


a x • F xM 


• • f xM 

•FxmJ 


( 31 ) 



( 32 ) 



G x = 


nv,™ H .r e , Ver ’ m0St °? h ° gonal P°'ynoniials are only orthogonal over a continuous range, not over discrete points To 

4. GRAM-SCHMIDT ORTH OG ONA LIZ ATI ON 
4.1 Summary of the Gram-Schmidt Orthogonalization process 

nnhocn? 6 Gram f Ch f midl or, hogonalization process takes M arbitrary vectors in an M dimensional space and creates M 
ve™ F,tl' : ™ * te ' S r °' ' hal SPa “ The G ' a "" Srf "” id ' Or.hogona)iMlion process coo * e“p“ ' £ “ 


j - i 3 k 

1 * ] ^XS ^XS 


(33) 


By letting D js represent the coefficient of 1 K (See Equation 34), then the onhogonalizalion of any vector F xj can 
be expressed as shown in Equation 35. y ’ J 1 


D = 

XJS ^ 

<X> • <J> 

xs xs 


j “ 1 

= **}- 

s = I 


(34) 


(35) 


4.2 Applying the Gram-Schmidt process to the vector formulation 

For the fitting process, e polynomial vectors, F x , through F,m (See Equation 22), are orthogonalized and replaced 
by the new orthogonal vectors, xi through ^ xM . Because the original fitting polynomials are replaced by orthogonalized 
polynomials, a different set of coefficients is found. These pew coefficients are named H, for j in 1 to M. Therefore, the 
original coefficients vector, G x , is replaced by the vector F x . The measured values can now be fitted in terms of the new 
vectors. 

“* = [*.. - (36) 

The system is then solved, as shown previously (See Equations 31 and 32), by multiplying both sides of the minimized least 
squares fit equation by the inverse of the resultant diagonal matrix. The T x coefficients are now known and can be used to 
represent the original surface. (See Equation 37) The process is repeated using the y-slope data to find r y 



( 37 ) 



5. FINDING 6* IN TERMS OF f x 

The polynomials fitted to the measured values are typically a standard set, such as the Zemike polynomials, 
erefore since the polynomials do not change, the coefficients found in the fitting process of the original polynomials (See 
Equation 32) can be usejl to describe the test surface and serve as a basis of comparison between other test surfaces. However, 
when the *i through xm were onhogonalized, the polynomials F x] through F xM were modified, creating a new set of 
polynomials, <D X , through <b xM , that are dependent upon the measured values. Thus, the T x coefficients that compose F, 
cannot serve as a basis of comparison, since evgry set_pf r x coefficients is related to a different set of polynomials. Because of 
this, it is necessary to find a means to convert to G Xf and F y to G y . 

_ • K' S, f P m findmg G,t m Ierms of r * is 10 express the orthogonalization process using matrix notation From 

Equation 35, the orthogonalization of all M vectors can be shown as follows. 
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D.m, d 
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5*.mI 


(38) 


The matrix of D x j, coefficients is called D x for simplicity. Solving Equation 38 for the orthogonal vectors, Equation 39 is 
found, where 1 is an identity matrix of size M*M. 


”_»T* 

1 

= d + D x )-i 

ft" 

1 

->T 



KmJ 


l^xM 


(39) 


Recalling Equation 23 and Equation 36, it is seen that two equivalent representations of the surface exist; one using the 

original polynomial vectors, and the other using the onhogonalized polynomial vectors. Therefore, this statement can be 
made: 


[ft - F*m] 


(40) 

Equation 40 is then rewritten using transposes (See Equation 41), apd the transpose pf both sides jaf the resultant equation is 
taken. This produces an expression for the matrix of through ^xm in terms of F xi through F lM (See Equation 42) 

K. - : 

= (i+d .)"[ 6, ...y 

(41) 

k, ^xm] ~ 

[f„ ...rJia+D.i-'jF 

(42) 

Now, Equation 40 is expressed using Equation 42. 



[ft Fxm] G * = 

>., ... F,J 1 c + D,)-'i T r, 

(43) 



r°' h A S l d ^n °! EqUa,i0n 43 are divided the matrix of po'ynomia] vectors, and G x i^fbund in terms of the D xjs coefficients and 
i . A e ain, this process is repeated using the y-slope data to find G y j n terms of r y 1 


G, = [(I + D,)-'] T f, 

6. WAVEFRONT RECONSTRUCTION FROM WAVEFRONT SLOPE 


( 44 ) 


As was stated in section 2, (INTERFEROGRAM ANALYSIS), a representation of the actual wavefront is obtained 
i ° 8 ° n s ’°P e representations. Therefore, it is necessary to solve for the weighting coefficients for the 

actual wavefront. This is done by using a linear combination of the orthogonal slope coefficients. (See Equation 45) The x- 

nT and y' 510 ^ c ° effi ^' en,s ^ multiplied by the weighting coefficients, H x and H respectively, which are of size M by M 
It should be noted that the scaling matrices are a function of the selected fitting polynomial set. 


G = H,G X + H y G y 


( 45 ) 


a n- ^ each e'ement of G depends only upon its corresponding slope elements (i.e. the polynomial set is orthogonal under 

' ,<>n •’ - en * if SC f In ° ^x anc * w *h he diagonal. However, if the fitting polynomial set is not orthogonal 

under differentiation, then fitting errors will occur, the Zemike polynomials are not orthogonal under differentiation; they 
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Term 

F(x,y) 

F x (x,y) 

F y (x ; y) 

0 

1 

0 

0 

1 

X 

1 

0 

2 

y 

0 

1 

3 

2 ( x 2 +y 2 ) - 1 

4x 

4y 

4 

| x 2 -y 2 

2x 

-2y 

5 

2xy 

2y 

-2x 

6 

3(x 3 +xy 2 ) - 2x 

3 ( 3x 2 +y 2 ) - 2 

6xy 

7 

3 {x 2 y+y 3 } - 2y 

6xy 

3 fx 2 + 3y 2 ) - 2 

8 

6 { x 2 +y 2 } 2 - 6 ( x 2 +y 2 ) 2 + 1 

24 (x 3 +xy 2 ) - 12x 

24 (x 2 y+y 3 ) - 12y 


There are two solutions to the problem of cross-correlation. The first is to orthogonalize the fitting polynomial sets, F 
and F y . The second option is to use a different basis set which is orthogonal under differentiation. The standard Zemike 
polynomials are such a basis set; they can be linearly combined to form their own derivatives. This is illustrated for the first 
nine Zernike polynomials: 




Term 

F(x,y) 

F x ( x , y ) 

F y (x,y) 

Zo 

1 

NA 

NA 

Zi 

X 

Z 0 

NA 

Z 2 

y 

NA 

Z 0 

Zj 

2 ( x 2 +y 2 ) - 1 


4Z 2 

Z 4 

x 2 - y 2 

2Zj 

-2Zj 

Z 5 

2xy 

2Z 2 

2Z X 

Z 6 

3 ( x 3 +xy 2 ) - 2x 

3 ( Z 3 +Z 4 ) + Z 0 

3Z 5 

Z 7 

3(x 2 y+y 3 ) - 2y 

3Z 5 

3(Z 3 -Z 4 ) ♦ Z 0 

Z B 

6(x 2 +y 2 ) 2 - 6(x 2 +y 2 ) 2 + 1 

8Z € + 4Zj 

8 Z 7 + 4 Z j 


Table 2 : The first nine Zernikes and combinations needed to create the x- and y-derivatives . 

In a similar manner, the weighting coefficients for a wavefront can be constructed from the weighting coefficients of 
the slope representation. This is illustrated for the first nine terms in Table 3. 


Term 

G 

Gx 

G y 

Go 

NA 

G x0 

G y0 

Gi 

G x0 " ( G x4 +G x3 +G y5> 

Gxl 

G yi 

g 2 

G y0 - (G x c>+Gy 4 + Gy3 ) /9 

G 2 

G y2 

G 3 

f G xl‘ fG y2)/® “ (G x g+Gy7)/16 

G x3 

G y 3 

g 4 

(G xl +G y2 )/4 - (G x6 -G y7 )/16 

g x4 

G y4 

g 5 

(G x2 +G y i ) /4 - (G x7 -G y6 ) /16 

G x5 

G yS 

g 6 

(G x4 + G x3+ G y5 )/9 - G xg /15 

G x6 

G V6 

g 7 

{ G x 5 “Gy4 _f Gy3 ) / 9 - Gyg/15 

G x7 

G y 7 

G 8 

(G x6 +G y7 ) /16 

G x8 

CD 

>* 
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Table 3: Weighting coefficients for the wavefront can be created from the slope coefficients. 
Recalling Equation 45, the information in Table 3 can be reduced to the scaling matrices, H x and H y as shown: 
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( 46 ) 


(47) 


7. CONCLUSION 

Ronchi shear interferometry is an optical testing technique similar to shearing interferometry. Information about the 
slope of a test surface is extracted from interferograms. This is done with three main steps. First, the slope data for two 
orthogonal directions is sampled, usually by one of two main methods: digitization or phase-measuring interferometry. 
Second, the sampled data points are ordered. And third, the ordered data is fitted to a the derivatives of a polynomial set, via a 
least squares method and Gram-Schmidt orthogonalization. A linear combination of the slope coefficients is then found to 
construct the wavefront from the slope representation. 
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ABSTRACT 

Interferometry is an optical testing technique based on the interference of light. Fringes are formed when the 
Optical Path Difference (OPD) between a reference beam and an object beam is an integral multiple of the illuminating 
wavelength. This OPD is extracted through the process of sampling, ordering and interpolating. This paper develops a 
linear-algebra vector notation model of the interferogram sampling and interpolation process. 

1. INTRODUCTION 

Whenever two wavefronts exist at the same point in space, they superimpose, to create an interference pattern. 
The irradiance at any point, E(x,y), as determined by the Optical Path Difference (OPD) between the two wavefronts is 
described by: 


E(x,y ) = £ 0 [1 + COS (JfcOPD) ] 


where 


2n 


* " T 

The two dimensional irradiance pattern created by the optical path difference at all points is called an irradiance pattern. 

The challenge is to extract a representation of the test surface from such an interferogram. This is accomplished 
by a three step process: data sampling, ordering, and interpolating. 


2. INTERFEROGRAM ANALYSIS 

Because an interferogram has nearly infinite resolution, the actual number of points that can be measured is 
limited only by the sampling technique. There are two main methods for acquiring the data: digitalization, and phase- 
measuring interferometry. In digitalization, data points are acquired by sampling along the center of either the light or 
dark fringes. Usually, the interferogram is photographed and the image is placed onto a digitizing tablet where the data 
points are digitalized by a human operator. Alternatively, the image can be acquired with a video camera and digitalized 
using image-processing techniques. Phase-measuring interferometry, by contrast, is an electronic process where a video 
camera samples the wavefront phase at every point within its detector array - thus creating a high density, uniform grid of 
data. An illustration to show the difference in the two techniques is shown in Figure 1 and Figure 2: 
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Fig. l.An illustration of tbe placement of data points 
when an interferogram is sampled by digitalization 



Fig. 2. An illustration of tbe placement of data points when an 
interferogram is sampled by Phase-Measuring Interferometry 

For this paper, it is assumed that the data is sampled by manual digitalization along the dark fringes. The data is 
recorded in the form (m,x,y) where (x,y) gives the position of the point and (m) is the order of the sampled fringe - which, 
as will be discussed in the next paragraph, is related to the height of the surface in units of 2n. 

Once the data has been sampled, it must be ordered. Because the irradiance pattern is proportional to 
cos(k*OPD) the data is in units of modulo 2tt or modulo X. Thus, even though a point in the interferogram represents a 
quantity greater than 2rt, no measured quantity will be larger than 2tt. For instance, an (x,y) point may have an (m) that is 
of magnitude 9n. The measured (m) will equal 0.5rc, which is found by taking the magnitude modulo 2n. Therefore, the 
magnitude of the data must be artificially reconstructed by ordering the data. This is done by assigning integer order 
number values to each series of measured data points. For example, all data points along a given fringe will all have the 
same order and each adjacent fringe has an incrementally larger order number. 

For manually digitalized data, this ordering is accomplished by the trained operator. For video digitalization, this 
is accomplished by image processing techniques. And for phase-measuring, this ordering is accomplished by a process 
called phase unwrapping. Once the data has been sampled and ordered, it must be interpolated to a polynomial set to 
represent the shape of the test surface, generally, this is accomplished using a least squares method with a Gram-Schmidt 
orthogonalization. From this fit several questions can be answered: how flat is the surface, are the features bumps or 
holes, are the features concave or convex, and what aberrations are present 

Traditionally, Zemike polynomials have been used for this representation for several reasons. They look like 
classical Seidel Aberrations. They provide a good fit to the circular aperture interferogram of most optical components. 
And, because they are orthogonal, they are easy to manipulate. It is possible to add or subtract coefficients without 
seriously affecting the quality of tbe fit 

To help visualize this three step process, two graphs are shown in Figure 3 and Figure 4. Figure 3 shows a one 
dimensional parabolic wavefront and the irradiance pattern due to the interference of the parabolic wave with a planar 
reference wave. This represents a slice of what would be a two dimensional interferogram. The pluses on the irradiance 
pattern represent the sampled points on the dark fringes. The numbers below the sampled points are the integer ordering. 
Figure 4 is a plot of the ordered points. It is easily seen that a parabolic function could be interpolated to those points to 
represent the test wave. 




Fig. 3. AID slice of a parabolic wavefront (top) and the irradiance pattern 
formed when this slice interferes with a plane wavefront (bottom). 

The pluses represent data points sampled at iiradiance minima. 
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3. VECTOR FORMULATION 

The purpose of this section is to describe simply and clearly the analysis process of an interferogram using 
general linear-algebra vector notation, and to develop a generalized mathematical model with which to explore the effects 
of: sampling parameters, fitting polynomials, and the interpolation algorithm on the accuracy of the reconstruction of the 
test surface. The remainder of this paper presents a general vector based formulation for the interpolation process. 



Any point in the interferogram, Z(x,y), can be represented by a series of M polynomials, F(x,y), and their 
coefficients, G. 


M 


z (*, y) * X °j F j ^ 
; = i 


Alternatively, the single point, Z(x,y), can also be expressed as a matrix of polynomials multiplied by a matrix of 
coefficients. 


[z(x,y)] = [fjUy) ••• f A/ (*’?)] 



Thus, if a sampled interferogram can be thought of as a series of N measured data points. 



v y l ),Z(x 2 , y 2 ) 



Then it follows that the complete set of all N data points can be represented as a polynomial matrix multiplied by the 
coefficients matrix. 
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In order to simplify this notation, the data points can be considered to be a vector of length N, a given 
polynomial evaluated at all N data points is a vector of length N, and the coefficients are a vector of length M. 
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Using this notation, the polynomial matrix can be expressed using F vectors. 

3 * [*1 ... t u | ( 3 ) 

Thus, rewriting equation 1, the surface can be represented as: 

i « 3G 

The polyno mials can be fit using a least squares method, where S is defined as the sum of the square of the 
difference between the data points and the fitted polynomials. If a perfect fit were possible, then S would be equal to zero, 
because there would be no difference between the data points and the representing polynomials evaluated at those points. 
However, there are always differences between the real surface and its representation. Therefore, S will be non-zero. The 
function of a least squares fit is to find the coefficients, for a given set of polynomials, which minimizes S. This process 
can be described using vector notation 

J »|[ 4 - 33]| 2 


= minimum 
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or as more commonly expressed using summations: 
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The polynomial coefficients are found by taking the derivative of S with respect to each coefficient and setting 
the result equal to zero. For example, minimizing S for G k gives: 

N N M 

X Z (V F k < V y r> = X X C /y ( V F * ( V 

r = 1 r= 1;= 1 


The result of minimizing S for all coefficients can be expressed using summations in matrix form as: 
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This expression can be greatly simplified with vector notation to give: 
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where the dot product of the vectors Fj and Fj gives the sum of the products of both polynomials over every data point. 

N 

p i • p j = X F i (v y r> F j y J 

r = 1 

Likewise the dot product of the vectors Z and Fj gives the sum of the products evaluated at all data points. 

N 

Z*Fj = X Z ( x r> y r> F j ^ V y J 

r = 1 

Finally, the minimization of S can also be written using the polynomial matrix and its transpose: 

3 T i = 3 r 3S 

where the transpose of a matrix is found by rotating its columns so that they become its rows. In the case of the 
polynomial matrix (See Equ. 3), the columns are composed of polynomial vectors. The polynomials vectors are rotated to 
become the rows of the polynomial matrix, as shown below: 



S T = 


Fm\ 


If the F vectors are orthogonal then the minimization process produces a diagonal matrix 
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and the coefficients can be found by multiplying both sides of the equation by the transpose of the diagonal matrix. 


G = 



( 5 ) 


However, most orthogonal polynomials are only orthogonal over a continuous range, not over discrete points. 
Thus, the solution process is not as easy as it first seems. To overcome this limitation, There are three possible 
approaches. The first is to assume that the polynomials are orthogonal. This may work if there is enough sampled data to 
approximate a continuous range, but it can introduce error. The second is to take the transpose of the non-diagonal matrix, 
but this can be numerically difficult and inaccurate. And third, the polynomials can be orthogonalized. Typically, this is 
accomplished using a Gram-Schmidt technique. Although this introduces extra work, it should pay for itself in the ease 
and accuracy it returns in the final solution. 

The Gram-Schmidt orthogonalization process takes N arbitrary vectors in an N dimensional space (N-space) and 
creates N ortho-normal vectors that form a basis for that space. Ortho-normal vectors are vectors which are orthogonal to 
each other, and are normalized, or have a unit length. A basis for an N-space is essentially a coordinate system for that 
space, since the basis vectors are linearly independent and all other vectors in that space may be formed from linear 
combinations of the basis vectors. 

Consider two arbitrary vectors, Fj and F 2 , which are wanted to be replaced by their orthogonal counterparts, <t>i 

and <X> 2 . 



To begin the process, the first vector is normalized by dividing by its magnitude.The result is defined as its orthogonalized 
counterpart. 
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After the first vector is orthogonalized, all other vectors are orthogonalized, in a systematic manner, by projecting the 
current vector onto each preceding orthogonalized vector and subtracting that projection from the vector leaving the 
orthogonal component 



This component is then normalized by dividing by its magnitude, and renaming it to its orthogonal counterpart. This 

’ .. t’o 

<J>/» s 


\F’r 



1 


process is repeated for each polynomial vector in the system. 

For the interpolation process, the polynomial vectors, F (See Equ. 2), are orthogonalized and replaced by the 
new vectors, <J>. Because the original fitting polynomials, F, are replaced by <I>, a different set of coefficients are found. 
Therefore, the original coefficients, G, are replaced by V. The system can now be rewritten using the new vectors, and 
because the vectors are orthogonal, the matrix is diagonal. 
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The system is then solved as shown previously by multiplying both sides by the transpose of the matrix. (See 
Equ. 4, 5) The T coefficients are now known and can be used to represent the original surface. 







4 . CONCLUSION 


Interferometry is an optical testing technique in which information about a test surface is extracted from an 
interferogram. This is done with three main steps. First, the data is sampled, usually by one of two main methods: 
digitalization or phase-measuring interferometry. Second, the sampled data points are ordered. And third, the ordered 
data is interpolated to a polynomial set. Typically, Zemike polynomials are used and are fitted with a least squares method 
and Gram-Schmidt orthogonalization. 

In order to describe the interpolating process, a generalized linear-algebra vector notation has been developed in 
which the data points, fitting polynomials, and coefficients are viewed as vectors. A least squares fit is applied to the 
vector system representing the test surface. The polynomial vectors are ortho-normalized via the Gram-Schmidt process 
and then system is solved for the coefficients vector. 
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Surface Deformation Measuring Instrument 

for the 

Surface Tension Driven Convection Experiment (STDCE-2) 


by: 

H. Philip Stahl, Ph.D. 



STDCE 


The Surface Tension Driven Convection Experiment 
(STDCE) is a fundamental fluid physics experiment 
designed to provide quantitative data on the 
thermocapillary flow of fluid under the influence of an 
increased localized surface temperature. 

STDCE flew on the Space Shuttle Columbia in the First 
United States Microgravity Laboratory (USML-1) in June 
1992. 

The second flight of this experiment (STDCE-2) is 
scheduled for flight in the Fall of 1995. 



SPACE EXPERIMENT CONFIGURATIONS LABORATORY APPARATUS SCHEMATIC SPACE EXPERIMENT APPARATUS 
NUMERICAL DEVELOPMENT AND EXPERIMENTAL VERIFICATION OF THERMOCAPILLARY FLOW IN REDUCED GRAVITY 
Fig. 1 Surface Tension Driven Convection Experiment. 





Motivation 


On Earth, buoyancy flow caused by gravity limits our 
ability to: 

grow semi-conductor crystals, 
grow biological crystals, and 
solidify metal alloys. 

by introducing imperfections into the lattice structure. 

Also, gravity causes lattice structures to ’sag 5 and limits 
the size of crystal which can form in suspension. 

In microgravity, buoyancy flow is reduced to 10' 6 , 
offering great promise for material processing in space. 

However, thermocapillary flow still exists and must be 
understood. 



Floating Zone Crystal Growth 
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FLOATING-ZONE SIMULATION 




FLOATING-ZONE SIMULATION 



TEMPERATURE OSCILLATION PATTERN IN LIOUIO 



Surface Tension 


Surface tension is a property of a liquid which tends to 
minimize its surface area. For, example, surface tension 
causes nearly spherical drops to form when rain falls. 

On Earth, gravity forces usually overpower surface 
tension forces, and liquids form flat surfaces. 


Convection 


Convection refers to the flow of a fluid. Natural 
convection is driven by gravity: colder, heavier fluid is 
pulled down, causing warmer, lighter fluid to rise. 

On Earth, this buoyancy driven convection occurs 
whenever there is a temperature difference within a fluid. 

Convection can also be driven by surface tension. As its 
temperature is increased, a liquid has a lower surface 
tension. Warmer surface liquid is pulled toward colder 
surface liquid, which has a higher surface tension. As the 
warmer liquid moves across the surface, it draws more 
liquid up to the surface, creating a convection current. 

This surface tension driven convection, also called 
thermocapillary flow, occurs whenever there is a 
temperature difference across the surface of a liquid. 


BASIC FLOW STRUCTURE 
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driven by thermocapillarity and 
viscosity 

driven by pressure gradient 



Thermocapillary Flow 


Thermocapillary flow has two states: 

Steady State, and 
Oscillatory. 

Steady flows are two-dimensional. They are the same for 
any symmetrical observer. 

Oscillatory flows are three-dimensional. 

Surface deformations produced by steady state 
thermocapillary flow may produce oscillatory flow. 

Oscillatory flow may be detrimental to the crystal growth 
process and needs to be studied. 


TIME LAG IN RETURN FLOW RESPONSE 



TOP VIEW 
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Return flow Surface flow 



CIRCULAR CONTAINER . 




Science Objectives 


The objective of STDCE is to determine the extent and 
nature of thermocapillary flow as a function of: 

heating mode and level, and 
liquid free-surface shape. 

The objective of STDCE-2 is to determine the onset 
conditions for and nature of oscillatory flow as a function 
of: 


heating mode and level, 
liquid free-surface shape, and 
container size. 

To accomplish these objectives, three parameters must be 
monitored: 

surface temperature distribution, 
bulk fluid flow, and 
surface deformation. 







STDCE Configuration 

The STDCE-2 hardware consists of five basic modules: 
Fluid Chamber, 

Heat Source (C0 2 Laser or Cartridge) 

Infrared Imager, 

Flow Visualization, and 
Surface Deformation. 

In STDCE, the fluid was in a single 10 cm diameter 
container. In STDCE-2, there will be three containers: 

1.2, 2, and 3 cm. 

The fluid is heated by illuminating the surface with a C0 2 
laser beam or a cartridge heater. 

Surface temperature is measured by a HgCdTe thermal 
imager. 

Bulk fluid flow is measured by monitoring the motion of 
microscopic particles which have been uniformly mixing 
into the fluid. This motion is tracked by shinning a 
’sheet’ of light through the fluid and observing the light 
scattered from the particles. 

STDCE had no means for measuring surface deformation. 
But, STDCE-2 will. 
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Figure 21.— Meniscus at edge of test cell. 
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Figure 22.— Comparison of numerical and IR imager surface temperature 
distribution for constant flux experiment with 0.90 , p = 2.1 W, 
T =22.7 *C, and beam diameter = 5 mm. Numerical calculation used 
Ma= I.27X 10 5 and Ra=4.87xl0 7 . 



Figure 23.— Constant flux experiment IR image for T = 38.7 ’C, R - 50, < — 0.90, p = 2.7 W, and beam diameter = 5 mm. 
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Figure 3 Schematic view of the light sheet illumination, fluid flow reservoir, and data 
acquisition and processing system. The tubular heater in the reservoir drives 
the flow. The frame-grabber board digitizes a sequence of video fields 
separated by AT, which is defineddn integral multiples of 1/60 second. 
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Combined data sequences #1 and #2. A large dynamic range of velocities 
have been recorded. Many sets of tail to head velocity vectors are observed, 
indicating particles which were tracked over the entire sampling period. 
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Interpolated 2-D -velocity vector plot. The combined data from sequences 
#1 and #2 have been interpolated over a 32*32 grid. The interpolated data 
clearly show the counter rotating flow pattern. 









Surface Deformation 


On Earth it has been observed that both steady-state and 
oscillatory thermocapillary flow states produce 
deformations of the air/oil free surface. 

Steady state flow produces two characteristic steady state 
surface deformations: 

a ’Pimple’ (bump), and 
a ’Dimple’ (hole). 

Pimples are produced when the surface is exposed to an 
external temperature source over a small area. 

Dimples are produced when the same temperature input is 
spread out over a larger surface area. 

Oscillatory flow produces both standing wave and rotary 
surface deformations. 
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Purpose 


The purpose of this project is to design, build, calibrate, 
and fully characterize one or more non-contact (i.e. 
optical) prototype instruments for measuring the surface 
deformation of an air-oil free-surface under the influence 
of a localized thermal load (i.e. C02 laser beam or 
electric heater) which can be integrated into the 
anticipated re-flight of the Surface Tension Driven 
Convection Experiment (STDCE-2) 

While initially these measurements will be performed in a 
laboratory environment, the instrument must be capable 
of operating in a Space Transport System (STS) 
environment for the United States Microgravity 
Laboratory flights USML-2 and/or USML-3. 



Measurement Approach 


Because the surface is a fluid, it cannot be contoured 
mechanically. It must be contoured optically. 

There are two ways to optically contour a large 
two-dimensional reflective surface: 

interferometrically or 
geometrically. 

Generally, it is easier and more accurate to use an 
interferometer than to perform a geometrical test because 
interferometers measure a surface’s shape while most 
geometrical tests measure a surface’s slope. 

Unfortunately, because of the large surface deformations, 
conventional visible interferometry is not practical. The 
fringes will be too dense for the camera to resolve. 

One solution to this problem is to use an infrared 
interferometer, but such a technique might interfere with 
the localized heating process. 

Therefore, classical geometrical techniques were 
evaluated: Foucault knife-edge test, wire test, Hartmann 
screen test, Ronchi test, grating interferometry, or 
reflective moire. 



Geometrical Techniques 


Geometrical techniques measure the surface’s slope 
instead of its height. Thus, there are specific advantages 
and disadvantages: 

Advantage: they are insensitive to tilting of the 
free-surface cause by G-Jitter. 

Disadvantage: two orthogonal data sets are required to 
reconstruct its shape. These two data sets can be either 
x-slope and y-slope or r-slope and 0-slope. 



Surface Symmetry 


A fundamental question which has a profound impact 
upon the selection of the final approach is whether or not 
the surface is rotationally symmetric. 

If the surface is rotationally symmetric, then it may be 
possible to satisfy the mission’s science requirements with 
just x-slope or r-slope. 

But, if the surface is non-rotationally symmetric, then it is 
necessary to acquire two orthogonal data sets. 

The surface deformations are not rotationally symmetric. 

But, their outer edge is at a constant height. 



Candidate Techniques 


Three techniques which provide data over the entire 
aperture were evaluated: 


Ronchi grating, 
Hartmann test, and 
Projected grid. 



Projected Grid 


In the projected grid approach, a linear or crossed grid 
pattern is placed directly in front of an extended light 
source (illuminated diffuser screen). 

This pattern is then imaged onto the reflective sample 
surface. 

The reflected wavefront is re-imaged by another optical 
system onto the video camera. 

The projected grid approach was base lined early but 
eliminated because it lacked sensitivity, dynamic range, 
and flexible spatial sampling. 



Ronchi Grating 


In the Ronchi grating test, a collimated beam of light is 
reflected off of the sample surface. 

The reflected wavefront is focused through a grating to 
form a Ronchigram which is imaged onto a video camera. 

The Ronchi approach has many advantages. By selection 
of grating it has variable sensitivity and dynamic range. 
By location of grating it has flexible spatial sampling. 

Unfortunately, it does not provide simultaneous 
orthogonal slope data sets. 




Observation 




Hartmann Screen Test 


The Hartmann screen test can provide both x- and y -slope 
information simultaneously. 

Unfortunately, it has limited spatial resolution. 

And, for this application, it is not practical because the 
screen would obscure the upper surface flow visualization 
function. 

A practical implementation of the Hartmann test is to use 
crossed gratings in a Ronchi test. However, the complex 
patterns produced are difficult to analyze. 



Two Channel Ronchi 


Another way to obtain simultaneous x- and y-slope data 
is with a two channel Ronchi approach. 

For the two-channel Ronchi approach, a collimated beam 
of light (laser or white) is projected onto the reflective 
sample surface. 

The reflected beam is collected, split into separate x- and 
y-slope measurement channels, passed through vertical 
and horizontal gratings, and projected onto an observation 
screen for viewing by a video camera. 

A two channel Ronchi instrument was demonstrated 
which simultaneously acquires x- and y-slope data for an 
input wavefront. Unfortunately, because of packaging, 
cost, and schedule issues, a two channel system could not 
be implemented. 





The Ronchi Test 


The Ronchi test is performed by observing the shadow 
pattern produced when a focused beam of light is passed 
through a periodic structure, such as a line grating. 

If the beam is ideal, the shadow pattern looks like the 
grating. But, if the beam is deformed, then the pattern is 
deformed. 

A Ronchigram is obtained when the shadow pattern is 
recorded at an image of the surface under test, such that 
it is superimposed (’painted’) upon the image. 

Information about the surface is obtained by correlating 
shadow lines with physical locations on the surface. The 
distance between adjacent shadows indicates how much 
the surface slope changes between those two locations. 

The amount of slope required to go from one shadow to 
its neighbor is called the equivalent wavelength. 






The Ronchi Test 


The Ronchi test can be explained as either a geometric 
ray or an interferometric test. 

From a geometric ray point of view, the fringes in a 
Ronchigram are the result of the deviation of a ray from 
its ideal path because of slope errors. 

From the interferometric point of view, the fringes in the 
Ronchigram are produced from the interference between 
overlapping diffraction orders. The overlapping 
wavefronts produce a shearing of the original wavefront. 

Therefore the Ronchi test can be modeled as a sheared 
wavefront interferometer. 



Advantages of the Ronchi Test 

Measures Slope: It is insensitive to G-jitter induced 
tilting of the oil surface. 

Self-Referencing: It is proportional to the difference 
between the test wavefront and itself displaced laterally. 

Common Path: The wavefront and its sheared self travel 
the same optical path, therefore, the test is insensitive to 
atmospheric turbulence and mechanical vibration. 

Known and Variable Sensitivity (Equivalent Wavelength): 
Measurement sensitivity is determined wholly by 
geometric factors (grating period and beam F/#) which 
can be varied in a controlled manner. 

Known and Variable Spatial Sampling: The number of 
spatial sample points across the aperture can be controlled 
by displacing the grating a known distance from focus. 



Disadvantages of the Ronchi Test 

Slope patterns are difficult to interpret visually. 

Two orthogonal data sets (x-, y-slope; or r-, 0-slope) are 
required to completely characterize a given wavefront. 

Orthogonal slope information can be obtained with a 
single-channel system using a crossed-grating, but the 
patterns are even more difficult to interpret and analyze. 


Shear Interferometry 


Whenever two coherent wavefronts are at the same point 
in space at the same time, they superimpose to form an 
interference pattern called an interferogram. 


The irradiance at any point, E(x,y), cause by this 
superposition can be expressed as a function of 
measurement wavelength and optical path difference 
(OPD) between the two interfering wavefronts. 


E(xy) 



1 + cos 


f 2n OPD > 

V A, ; 


\ 
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where : X = measurement wavelength 


For a shear interferogram, the OPD is given by the 
difference between the wavefront and its sheared self. 

OPD = W(x>y) - W(x + Ax j) = A W(xy) 
where : Ax = shear distance 

A W{xy) = differential wavefront 



Shear Interferogram 

• The irradiance for a sheared interferogram is given by: 

• Relationship between differential wavefront and slope is: 

AWfoy) = 5 Ax = a(xy) Lx 

bx 

where : a(xyO = wavefront slope 

• Inserting this into the irradiance equation gives: 

/ 

E(x,y) = E 0 1 + cos 
1 

where : X = — = equivalent measurement wavelength 
q Ajc 




Ronchigram 


Similarly, a Ronchigram is described by: 


E(*y) = E , 


1 + COS 
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where: a' = 


bW 

bx' 


= Normalized Slope 


x' = - = Normalized Pupil Coordinate 
r 


r = Pupil Radius 


X = — 
eq 2 F/# 

d = Grating Spacing 

F/# = 7&sf Ream F/# 


Thus, with proper scaling, conventional interferogram 
analysis provides a direct measurement of the slope. 

Finally, if the grating is shifted forward of focus, the 
irradiance pattern becomes: 
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Ronchi Theory 



Screen 


• The wavefront in the front focal plane of the 
collecting lens is 


«(*i0’ 1 ) 


• The wavefront in the back focal plane of the 
collecting lens (without the grating) is the Fourier 
transform of the wavefront in the front focal plane. 





Ronchi Theory Continued 



• The wavefront is modulated by the transmission 
function of the grating. 


«( X 2 >/ 2 ) = u ^y/x? ^ 2 ) 


• At the observation screen, the wavefront is given by 

«(* 3 ,y 3 ) = ^{u(y 2 ,y 2 )} 


u(x 3 ,y 3 ) 


, 1) 


Xz Xz 


When z = f the wavefront at the observation screen 
is the original wavefront convolved with the Fourier 
transform of the grating transmission function. 




Sine Wave Grating 


• A sine wave grating is a series of transparent and 
opaque straight line bands of the following form. 



t(x 2 ,y 2 ) = 1 + cos(-) 

a 


• The Fourier transform of the sine wave function is 
three delta functions. 

+ cos(— )} = &(*) + 6(jc - d) + 6(x + d) 
d 



Sine Wave Grating 

y3 



E = E a + E b + E c + {u a u' b + u b u*) + (u a u* + u c u*) + (u b u* + u c u b ) 


The cross terms (overlapping regions) are the terms that 
interfere to produce the Ronchigram. 


Conceptual Desi 



Aeq = d / [2(f/#)] x' = x / hi 



Flight System 


Based upon managerial issues (packaging, cost, and 
schedule) and science issues (sensitivity, dynamic range, 
and spatial sampling) a single channel Ronchi system was 
selected. 

The managerial constraints and science requirements were 
then appropriately modified. 



Surface Deformation Science Requirements 


Observe the entire surface for three different sized 
chambers: 12 mm, 20 mm, and 30 mm. 

Measure slope as small as 5 jum/mm. 

Measure slope as great as 30 fimjmm without vignetting. 

Spatially sample the surface at least 10 times per 
diameter. 

Temporally sample the surface at least 30 times per 
second. 

Reconstruct the surface shape at least along a diameter. 
Indicate when the oil is filled to its flat surface position. 



Science Requirement 
Field of View 

The surface deformation measuring apparatus must be 
able to view the entire surface area of each fluid chamber 
at an optimum camera spatial resolution. 

The chamber sizes are: 

30 mm 
20 mm 
12 mm 

The apparatus’s imaging system should have at least three 
different fields of view which can be used to view any of 
the three chambers. 



Science Requirement 
Measurable Surface Features 

The apparatus must be able to measure three different types of 
surface features: 

• A central bump or ’pimple’, 

• A central hole or ’dimple’, and 

• A large-area temporal oscillatory motion. 
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Science Requirements 
Central Bump 


The central bump is assumed to be spherical in shape 
with a constant radius of curvature, given by: 


W(y) = A 


i - y- 


where: A = Maximum Vertical Height, and 

r = Maximum Lateral Radius. 

The maximum measurable surface slope of the central 
bump occurs at y = r and is given by: 


bW(y) 

6y 


2 A 


I max 


The maximum measurable surface curvature of the central 
bump is given by: 


6 2 W(y) 

by 2 


2 A 


max 




Science Requirements 
Centra] Hole 


The central hole is assumed to be Gaussian in shape and 
is given by: 


“71 


where: 


A 

r 


W(y) = - A e' r 1 

, Maximum Vertical Height, and 
Maximum Lateral Radius. 


The maximum measurable surface slope of the central 
hole occurs at y 2 = r 2 /2k and is given by: 


bW(y) 

6y 


y/2 71 A ( 2 ) 1.5 A 


1 max 


The maximum measurable surface curvature of the central 
hole occurs at y = 0 and is given by: 


b 2 W(y) 

by 2 


2k A 


6.3 A 


1 max 







Science Requirements 

Maximum Resolvable Features 

The largest feature that the apparatus must resolve, 
whether it is a bump or a hole, is assumed to have the 
following properties: 

Vertical Height A = 250 jum 
Lateral Radius r = 5 mm 

The maximum assumed measurable surface slope and 
curvature for a bump are: 

Slope 100 //m/mm 

Curvature 20 //m/mm 2 

The maximum assumed measurable surface slope and 

curvature for a hole are: 

Slope 75 jum/mm 

Curvature 63 //m/mm 2 



Interferogram Analysis 


A minimum of 3 pixels per fringe are required to 
extract the correct information about the test signal. 


Rise, Run, and Feature Sign 


A good rule of thumb is to have 4 pixels per fringe. 


The curvature corresponding to one fringe per 4 
pixels can be calculated by: 


P = 


4ar 


where P = curvature , k = equivalentwavelength, 
a = pixel width , and r = radius of test beam . 


For = 60.48 um, r = 10 mm, and a = 19.5 um the 
curvature is calculated to be 77.5 um/mm 2 

Since the calculated curvature is larger than the 
required curvature of 63 um/mm 2 , the required 
curvature can be resolved and each fringe covers 
more than four pixels. 

Since each cycle of the test function is larger than 
four pixel widths the Nyquist rate Is met (§, >2§ 0 ). 



Science Requirement 

Minimum Resolvable Feature 

The smallest feature that the apparatus must resolve, 
whether it is a bump or a hole, is assumed to have the 
following properties: 

Vertical Height A = 10 jum 
Lateral Radius r = 5 mm 

The minimum assumed measurable surface slope and 
curvature for a bump are: 

Slope 4.0 //m/mm 

Curvature 0.8 //m/mm 2 

The minimum assumed measurable surface slope and 
curvature for a hole are: 


Slope 3.0 //m/mm 

Curvature 2.5 //m/mm 2 



Ronchi Measurement Performance 


The performance of the Ronchi system is completely 
specified by its equivalent wavelength: 

_ d 

eq 2 (F/#) 

Since the F/# of the measurement beam is fixed by the 
optical system, everything depends upon grating spacing. 

This dependance is extremely important. It allows for the 
measurement capability of the system to be varied on 
orbit by simply inserting different gratings into the beam. 

If the surface deformation is large, insert a coarse grating. 

If it is small, use a fine grating. 



0.5 cvcles/m ni 1.0 cycles/mm 2.0 cycles/mm 

Figure 14 Comparison of Slope Detectable with Different Gratings 





Performance Parameters 


Accuracy defines the uncertainty of a given measurement. 

For the current sampling system, the measurement 
accuracy is approximately X eq /8. 

A measurement system can report numbers to an arbitrary 
precision, but they may not be accurate or repeatable. 

The minimum measurable slope of a given system is 
defined to be twice its measurement accuracy. 

For the current sampling system the minimum measurable 
slope is approximately X eq /4. 

The maximum amount of slope which the system can 
measure is determined by how many shadow lines 
(fringes) the video camera can clearly resolve in the 
Ronchigram. 

For the current sampling system, the Ronchigram will be 
limited to approximately 16 lines pairs for a maximum 
slope of 16 



Table 1. Performance Summary for F/2.25 System 


Grating 

Grating 


Accuracy 

Minimum 

Maximum 

Frequency 

Line Space 

^“eq 

[am/mm] 


Slope 

Slope 

[lp/mm] 

[mm/lp] 

[um/mm] 

[am/mm] 

[am/mm] 

0.5 

2.0 

444.44 

55.56 

111.11 

7111 

1.0 

1.0 

222.22 

27.78 

55.56 

3556 

2.0 

0.5 

111.11 

13.89 

27.78 

1778 

2.5 

0.4 

88.89 

11.11 

22.22 

1422 

5.0 

0.2 

44.44 

5.56 

11.11 

711 

10.0 

0.1 

22.22 

2.78 

5.56 

356 











































Science Requirements 
Temporal Resolution 


The third required surface feature which must be 
measured is a large-area temporal oscillation. 

It is assumed that the period of this oscillation is on the 
order of 2 to 5 seconds. Thus, a standard video camera 
with a 60 Hz field rate is wholly adequate as the data 
acquisition device. 

However, if the temporal phenomena were to have a 
period of less than about 0.2 sec, then a non-standard 
high-speed video camera may be necessary. 

It is assumed that there may be high-speed temporal gitter 
in addition to the lower-speed temporal oscillation. The 
effect of this gitter is to reduce the contrast of the pattern. 

To eliminate gitter the camera will need a high speed 
shutter and the source will need sufficient output power 
to expose the camera. 



Science Requirement 
Fill Level Indication 

Silicon oil has a very low surface tension and spreads 
easily on any surface. To prevent it from flowing out of 
the container, the rim of the chamber has a sharp edge 
that ’pins’ the oil in place and is ’barrier’ coated. 

As the container is filled in low gravity, the liquid will 
form a deep spherical curvature. 

When completely filled, the surface will be flat. 

If over-filled, the surface will form a spherical dome. 

By placing a grating at a focal point, direct visual 
evidence of the oil’s surface state is provided to the 
Mission Specialist. 

When the oil surface is not flat, there is a line pattern. 

As the oil approaches a flat surface, the pattern spreads 
out (has fewer lines) until all the lines disappear when the 
surface is flat. If the container is under or overfilled, the 
lines will reappear and increase in number. 



Concave Surface 


Flat Surface 


Convex Surface 


Figure 9. Fill Level Indicator 





Flight System Design 

The free surface deformation measuring instrument was 
designed in modules. 

The flight system consists of six modules: 
light source, 
projection/imaging, 
polarization, 
pupil relay, 
measurement, and 
camera. 

Each module overcomes various engineering and/or 
packaging challenges while allowing the total system to 
meet the science requirements. 

Each module is independently assembled and bench 
aligned before integration into the final system. 

And, each module prevents unnecessary exposure of the 
crew to light radiation. 



Source and Beam Expansion Module 



A 





Source Module 


The source module provides a collimated beam of light 
which uniformly illuminates the object under test. 

This is accomplished with a laser diode and a Galilean 
beam expanding telescope. 

A laser diode was selected because it is a very-bright 
collimated source capable of surviving the launch. 
(White-light and LED sources were also considered.) 

A commercially available Galilean beam expander was 
selected because it is compact and does not have a spatial 
filter which could become misaligned during launch. 

The beam diameter is defined by the projection/imaging 
module’s entrance pupil to be 22 mm in diameter. 

The beam is magnified (30X) and truncated to balance 
uniformity and power. 



Conceptual Design 




Fluid Properties 


The object under test is a free surface of 2 cSt silicon oil. 

It has a refractive index of 1.39 

And, an irradiance reflection coefficient of 2.7%. 

Because of this low reflection coefficient, the source must 
have at least 10 mW of power for the video camera to 
observe the Ronchi pattern. 

Also, ghost reflections from the optical components can 
obscure the oil surface reflection. 



Illumination Geometry 

Because the free-oil surface is a specular reflector, we 
treat it as a plane mirror. 

A collimated illumination beam is required to implement 
the Ronchi test. 

Two illumination geometries were considered: 

Off-Axis 

On-Axis 

The on-axis geometry was selected for several reasons: 

Able to get lens closer to fluid surface, thus able to 
observe the highest slope errors. 

No off-axis projection distortion or cosine scaling. 

Easy to package. 

The on-axis geometry requires the beam to travel in two 
directions through the projection/imaging module. This 
can cause ghost reflection problems. 



Conceptual Des 
Off-Axis Illumination 




Conceptual Design 
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Surface Projection Imaging Module 





Conceptual Design Off Axis Illumination (Distance to First Lens) 



f >= I/2sin(0)[Dcos(0) - Max.] 
if D = 40 mm Max. = 30 mm 

f >= 41 mm 


Conceptual Design On Axis Illumination 



f >= 30 mm 



Conceptual Design Off Axis Illumination (Chamber Size vs. Vignetting) 



Container 



Conceptual Design Illumination Geometry Pupil Mapping 





Major Axis = (Max. Size) / cos (0) 

When projected onto the camera the axis will switch 


Polarization Module 


The polarization module eliminates ghost reflections and 
efficiently uses the source’s available optical power. 

Since the projection/imaging module is used in both 
directions and because the oil’s reflection coefficient is 
small, ghost reflections from the illuminated optical 
surfaces can seriously obscure the oil reflection. 

To eliminate this problem, the polarization of the beam is 
manipulated such that the ghost reflections are vertically 
polarized and the oil reflection is horizontally polarized. 

This is accomplished with two polarizers and a quarter 
wave plate. The first polarizer defines the polarization of 
the ghost reflections. The quarter wave plate rotates the 
oil reflection polarization by 90 degrees. And, the second 
polarizer blocks all of the ghost reflections and passes the 
oil reflection. 

To maximize the available power, a polarization beam 
splitter is used. All of the source’s vertical light is 
transmitted into the projection module. And, all of the 
returning horizontal light is reflected into the relay 
module. A conventional 50/50 beam splitter could be 
used, but it would throw away 75% of the available light. 



Conceptual Design 
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Surface Projection Imaging Module 







Figure 6. Polarization Module 



Projection/Imaging Module 

The projection/imaging module has two functions: 

illumination 

imaging. 

Its illumination function magnifies or de-magnifies the 
light beam to illuminate the oil surface for each of the 
three chamber sizes (12 mm, 20 mm, 30 mm). 

Its imaging function has two requirements: 

form an image of the different sized surfaces under 
test at a fixed sized pupil location, and 

pass a reflected wavefront with at least 30 jum/mm of 
slope without vignetting. 

The illumination and imaging functions determine the 
magnification properties of this module. 

The vignetting requirement determines its F/#. 

Packaging issues determines its clear aperture. 



Projection/Imaging Module 
Illumination Function 

The illumination function is accomplished with two 
AFOCAL lens pairs, an outer pair and an inner pair. 

The outer pair forms a Keplerian telescope with unit 
magnification. 

The inner pair forms a Galilean telescope which is 
positioned about the outer pair’s internal focus and 
functions as a pseudo field lens. This lens provides 
magnification/de-magnification when it is flipped. 

To illuminate the 20 mm chamber, remove the flip lens 
such that the 22 mm source beam is relayed unchanged. 

To illuminate the larger or smaller chambers, insert the 
flip lens either forward or backward to magnify or de- 
magnify the beam. The final illumination beam size 
depends upon the flip lens magnification factor. 

The flip lens is not symmetric about the focal point and 
does not rotate about the focal point. 



Conceptual Design Imaging Module 




0.5X Magnification 





Conceptual Design Imaging Module (2X, 3X) 
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Conceptual Design AFOCAL Flip Lens 
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Projection/Imaging Module 
Imaging Function 

Imaging of the different sized chambers into a fixed size 
pupil is accomplished by running the beam backwards 
through the system. 

The Keplerian telescope, without the flip lens, relays an 
image of the surface in the middle sized chamber. 

Inserting the flip lens relays images of the larger or 
smaller sized chambers. 

For proper imaging the surface must be in the front focal 
plane of the lens closest to the oil. 

The resultant image will be in the back focal plane of the 
last lens of the projection/imaging module. This image is 
transferred to the camera module by the pupil relay and 
measurement modules - forming a Ronchigram. 



Projection/Imaging Module 
Vignetting 

The projection/imaging module (as well as the pupil relay 
and measurement modules) must pass a maximum slope 
of 30 /rni/mrn without vignetting. 

This places an F /# requirement on these components. 

The faster their F/# the more slope they can pass. 

Initially, the optical system was designed to fully collect a 
wavefront reflected from a 250 jum high, 10 mm diameter 
spherical deformation. 

The actual design was determined by packaging and 
schedule issues. 

Packaging factors limit their maximum diameter and 
minimum focal length. 

Schedule require them to be commercially available. 



Conceptual Design Reflected Wavefront from Bump 
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= (r) A 2 / (2 6) = (5 mm)'2 / (2(0.250 mm)) = 50 mm 



Vignetting Requirement 

The selected components are F/2.25 (90 mm focal length, 
40 mm diameter) achromatic doublets with approximately 
2 waves of spherical aberration at 633 nm. 

This amount of aberration is acceptable given the 
magnitude of the anticipated oil surface deformations. 

With proper calibration, this error can be removed from 
the measurement. 

Please note that the 30 jum/mm specification is not 
satisfied for the 30 mm diameter chamber. 
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Pupil Relay Module 

The pupil relay module is a unit magnification AFOCAL 
system. 

Its original function was simply to extend the beam path 
such that the measurement module is at a location where 
the mission specialist can insert the gratings. 

Additionally, it serves as an alternative measurement 
module allowing for gratings to be inserted into parts of 
the beam not accessible with the measurement module. 

This is required to properly measure hole deformations. 




Expansion of Relay 
Module in Figure 3. 





Measurement Module 


The measurement module allows gratings of different line 
spacings to be inserted into different locations of the 
focused oil reflection forming a Ronchigram. 

Since the lens F/#’s are known, changing the grating line 
spacing changes the equivalent wavelength. 

The equivalent wavelength depends only upon the grating 
line spacing and is independent of where the grating is 
placed in the beam. 

Grating placement in the beam determines the number of 
spatial sample points across the beam. 

To insure at least 10 sample points per diameter, a coarse 
grating must be placed at a greater distance from the 
measurement module focus than a fine grating. 

The distance from focus for a given grating depends 
entirely upon its line spacing: 

z = 10 d (FI#) = 22.5 d 

To reconstruct a two-dimensional surface contour, 
separate horizontal and vertical gratings, or crossed 
gratings can be inserted into the beam. 

Radial or circular zone gratings may be tried. 




Measurement Module 
Expansion of Measurement 
Module in Figure 3. 
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0.5 cvcles/mm 1.0 cycles/mm 2.0 cycles/mm 

Figure 14 Comparison of Slope Detectable with Different Gratings 




Table 2. Distance from Focus for Spatial Sample Points 


Grating Grating Z for 10 Z for 15 Z for 20 

Frequen Line Samples Samples Samples 

cy Space per per per 

[lp/mm] [mm/lp] Diameter Diameter Diameter 



5.0 0.2 4.5 mm 6.75 mm 9.0 mm 

10.0 0.1 2.25 mm 3.375 mm 2.5 mm 




































Camera Module 


The Ronchigram produced by the measurement module is 
imaged onto a diffuser plate where a real image is 
formed. 

This image is viewed by a video camera and recorded on 
video tape for subsequent data analysis. 

The physical distance from the diffuser screen to the 
camera is determined by the focal length of the camera 
lens, the camera sensor format, and the size of the image 
on the diffuser screen: 


d = 


j + j hmage 


sensor 


i 


camera lens 


For the flight system, the diffuser image radius is 10 mm, 
the sensor radius is 2.2 mm, and the camera lens focal 
length is 28 mm. Thus, the stand-off distance is 155 mm. 


The camera is a 1/2 inch format RS-170 video camera. 



Calibration 


If the optical system were aberration free, the 
Ronchigram would be a series of perfect straight lines. 

However, given the total number of positive optical 
components in the optical system, such a pattern is 
impossible 

By measuring the beam reflected from a ’perfect’ 
reference-flat insert, after assembly, these errors can be 
characterized and removed from all measured data. 



Certification 


The prototype was certified by measuring a known 
amount of defocus aberration introduced by translating 
the collimating lens. 



8 (F/#) 2 


where', e, = Longitudinal Aberration 


Defocus was used because it is: 
easy to introduce, 
linear with translation, and 
easy to analyze. 


Additionally, the prototype measured diamond-turned 
samples provided by NASA. 

Measurement accuracy is approximately 5%. 

The flight instrument will be certified by measuring 
known test plates. 



Calibration results 


For each grating, a linear relationship was observed 
between the amount of lens translation and the amount of 
measured defocus. 

Multiplying each defocus value by the appropriate 
equivalent wavelength produces a single line of data. 

The average percent error between the data and the best 
fit line is ± 5%. 
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Grating Alignment 

The system is designed for the insertion of 
interchangeable gratings at various locations in the beam. 

Thus, grating alignment errors can be a problem. 

The data analysis was evaluated for its sensitivity to error 
in the position of one grating relative to its orthogonal 
partner: 

Rotation about the Surface Normal 
Rotation parallel to the Line Structure 
Rotation perpendicular to the Line Structure 
Translation along Surface Normal 
Translation parallel to Line Structure 
Translation perpendicular to Line Structure 

In general the Ronchi test is relatively insensitive to small 
grating placement errors. 

Except for translation along the surface normal. 





Figure 16 Vectors Corresponding to the Slope in the Ronchigrams in 
Figure 15. 





Figure 17 Ronchigrams with Grating Misalignment 
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Figure 18 Vectors Corresponding to the Slope in the Ronchigrams 
Figure 17. 







Grating Rotation 

Axis Parallel to Grating Structure 



Defocus — Astigmatism Coma 














Grating Translation 

Along Surface Normal 



Defocus — •— Astigmatism -e- Spherical 




Post-Flight Data Analysis 

The recorded Ronchigrams are analyzed after the flight to 
reconstruct the surface shape along a given diameter - and 
if possible over the entire surface. 

This is accomplished using standard interferogram 
analysis techniques to obtain a polynomial representation 
of each slope pattern. 

Surface shape along a diameter is obtained by integrating 
across the slope data. 

With orthogonal data sets, surface shape is obtained by a 
linear combination of Zemike polynomials. 



Conclusions 


A prototype Ronchi instrument to monitor: 

Steady State Deformation, 

Oscillatory Deformation, and 
Fill Level 

of a free-surface of oil has been developed, demonstrated, 
and characterized for STDCE-2. 

A single-channel Ronchi configuration has been selected 
as the flight instrument. 

The flight instrument is currently being designed 



Conclusions 


The Ronchi test is a shearing interferometer. 

The number of sheared wavefronts and the magnitude of 
the shear is dependant upon the characteristics of the 
grating. 

The fringe spacing (equivalent wavelength) in a 
Ronchigram is proportional to the spacing of the grating 
and the F/# of the test beam. 

The dynamic range and sensitivity of the Ronchi test can 
be changed by simply changing the grating spacing. 

The two channel Ronchi system has an accuracy of ± 5%. 

For a shearing system, the maximum measurable error is 
determined by the surface curvature. 
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